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NAPOMENA: Fotokopiranje ili umnozavanje na bilo koji nacin ili ponovno objavljivanje ove knjige - u celini 
ili u delovima - nije dozvoljeno bez prethodne izridite saglasnosti i pismenog odobrenja izdavaca. 


Predgovor 

Kao sto je poznato, pravac i brzina vetra u nekom mestu uslovljeni su 
raspodelom vazdusnog pritiska na teritoriji na kojoj se posmatrano mesto nalazi. 
Vazduh u horizontalnom pravcu struji uvek od oblasti visokog ka oblasti niskog 
vazdusnog pritiska, tj., od anticiklona ka ciklonu. Ovo strujanje se ne vrsi 
direktno u pravoj liniji iz razloga sto na jedan vazdusni delic, kod krivolinijskih 
izobara, deluju uglavnom cetiri sile koje imaju razlicite pravee, smerove i 
intenzitete. Te sile su: sila gradijenta vazdusnog pritiska, sila devijaeije, sila 
trenja i centrifugalna sila. Pod uticajem ovih sila delici vazduha se krecu po 
njihovoj rezultanti. Uglavnom vazdusne struje u prizemnim atmosferskim 
slojevima skrecu na severnoj polulopti udesno od pravea gradijenta pritiska, a na 
juznoj polulopti ulevo. 

Ako vazduh duze vremena miruje iznad neke oblasti on dobija izvesne 
fizicke osobine te oblasti u odnosu na temperaturu, vlaznost itd. Medutim, kada 
je vazduh, pod uticajem navedenih sila, primoran da se krece u horizontalnom 
pravcu, on te svoje osobine menja ili gubi i donosi oslabljene u oblasti u koje 
dolazi, odnosno preko kojih prelazi. Prema tome, horizontalna vazdusna 
strujanja iznad nekog mesta prouzrokuju neperiodske promene temperature, 
vlaznosti vazduha, oblacnosti itd., a u vezi sa ovim promenama izazivaju i 
izvesne meteoroloske pojave kao sto su kisa, sneg i dr. 

Kada se sve ovo ima u vidu, onda se moze, na osnovu proucavanja pravea i 
brzine vetra u nekom mestu ili oblasti, kao i proucavanjem zavisnosti ostalih 
meteoroloskih elemenata i pojava od pravea i brzine vetra, izvrsiti analiza 
vazdusnih masa koje prelaze preko tog mesta. 

Pod uticajem reljefa, narocito ako oblik i polozaj planina tome posebno 
pogoduju, moze doci do takvih poremecaja i modifikaeija opstih vazdusnih 
strujanja, koji mogu usloviti pojavu karakteristicnih lokalnih vetrova. Ovi 
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lokalni vetrovi javljaju se na zavetrenim stranama planina, na kojima vazduh 
prelazi u silazno kretanje. Orografske vetrove, s obzirom na njihovu izvesnu 
slicnost sa slapovima (vodopadima), nazivamo slapovitim vetrovima. Jedan od 
takvih vetrova je i kosava. 

Glavni uzroci kretanja vazdusnih masa u nasim krajevima, vezani su za 
promene u velikim akcionim centrima kao sto su Sibirski anticiklon, Islandska 
depresija, Azorski anticiklon i Mediteranska depresija. U zimskim mesecima 
dominantan uticaj na vreme kod nas imaju Sibirski anticiklon i Mediteranska 
depresija. Pri takvoj sinoptickoj situaciji izobare iznad naseg podrucja pruzaju se 
u smeru jug-sever. Prizemni vetar u toj situaciji je jugoistocni (kosava). Visinski 
vetar (koji duva duz izobara) duva iznad naseg podrucja kao juzni ili cak 
jugozapadni. U zimskim mesecima, kada dominira Sibirski anticiklon, na istoku 
je visok pritisak, koji blokira osnovno zapadno strujanje (koje uglavnom 
preovladuje u letnjim mesecima). Drugim recima, kosava je vezana za 
kontinentalni uticaj i uporedno prisustvo Mediteranske depresije. 

Kosava ima uticaja i na padavine, tj., nema pojave kise dok duva kosava. 
Tako kosava ima i izvestan prognosticki znacaj. 

U ovoj monografiji materija je grupisana u sedam glava. U uvodnoj glavi su 
prikazane dve sinopticke situacije koje uslovljavaju pojavu dva tipa kosave na 
teritorij i Srbije, kao i rezultati statisticke obrade vetra na teritoriji Beograda. 
Monografija obuhvata pored osnovnih karakteristika vetra kosave (Glava 2) i 
prostornu mezoanalizu dve kosavske situacije (Glava 3). Statisticka analiza 
maksimalnih brzina vetra kosave kao i udarnosti data je u Glavi 4. Promena 
brzine vetra, temperature i vlaznosti vazduha sa visinom u danima sa najduzom 
kosavom detaljno je opisana u Glavi 5. Spektralna analiza vetra kosave, koristeci 
Blekman-Takijev metod i metod sa vise prozora, kao i u novije vreme aktuelnu 
vejvlet transformaciju, prikazana je u Glavi 6. Razvoj prognostickih modela (na 
primer Eta modela) omogucio je i primenu ovih modela za simuliranje kosave, 
sto je detaljno prikazano u Glavi 7 za dve izabrane sinopticke situacije. 

Na kraju, drago nam je sto mozemo ovim putem da se zahvalimo prof. 
Mladenu Curicu i prof. Lazaru Lazicu koji su pazljivo procitali rukopis i dali 
korisne sugestije. 
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1 Uvod 


Posto skoro sve ljudske aktivnosti vise ili manje zavise od vremena, to su 
tokom zadnjih decenija dvadesetog veka meteorolozi najvise paznje posvetili 
upoznavanju i prognoziranju razlicitih klimatskih elemenata koji cine vreme. 
Tako, u meteorologiji mozemo govoriti o citavim sistemima koji se svojim 
nastankom, razvojem, kretanjem i iscezavanjem razlikuju od ostalog dela 
atmosfere. Po prostornim razmerama ove sisteme mozemo podeliti na 
mikrosisteme, mezosisteme i makrosisteme, kojima su se meteorolozi najvise 
bavili. 

Medutim, makrosistemi, koji su odgovorni za pojavu vetra kosave na vecem 
delu nase teritorije su veoma znacajni za razvoj vremena. 

Postojano strujanje vazduha iz oblasti veceg pritiska ka oblasti manjeg 
pritiska naziva se vetar. Vetar ima svoj pravac, smer i jacinu, pa je stoga 
vektorska velicina. Vetrovi mogu biti postojani kao sto su istocni i zapadni 
vetrovi ili povremeni kao sto su lokalni vetrovi. Lokalni vetrovi u nasim 
krajevima mogu biti vetar s kopna ili vetar s mora, burajugo i kosava, vetar koji 
je osnovna tema ove Monografije. 

Svi podaci o vetru (kosava) korisceni u ovoj Monografiji registrovani su u 
Beogradu na Meteoroloskoj Opservatoriji na Vracaru ((p=44°48' N, A,=20°28' E, 
H = 132 m) i Aeroloskoj Opservatoriji na Zelenom Brdu ((p=44°47' N, A,=20°32' 
E, H = 243 m), kao i na drugim meteoroloskim stanicama u kosavskom podrucju 
(juzni Banat, Podunavlje i Pomoravlje). 

Podaci o vetru se dobijaju osrednjavanjem sa anemografske trake. Primer 
umerene kosave zabelezen je na anemografskoj traci 20. novembra 1991. godine 
u Velikom Gradistu (slika 1). 
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Slika 1. Anemografska traka zabelezena u Velikom Gradistu 20. novembra 1991. godine. 


Statisticka obrada vetra uglavnom se svodi na analizu cestine pravaca vetra 
u prom i lima (%o), na osnovu koje se crta ruza vetra, kao i na analizu srednjih 
brzina vetra za svaki od 16 pravaca, na osnovu koje se crta histogram srednje 
brzine vetra tokom godine (slika 2). 
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Slika 2. Srednja cestina i srednja brzina vetra za 16 pravaca za stanicu Beograd-Opservatorija u 

periodu od 1888. do 1991. godine. 



U tabelama 1 i 2, koje su koriscene za crtanje slike 2, pokazane su destine 
16 pravaca vetra i tisina (%o) kao i srednje brzine vetra u Beogradu u periodu od 
1888. do 1991. godine. Pored znatnog broja tisina (183%o) za Beograd je 
karakteristicna pojava vetra iz jugoistocnog pravca-SE (145%o) koji je poznat 
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kao kosava (tabela 1). Najvecu cestinu ovaj vetar ima u hladnoj polovini godine 
(od 1 58%o do 189%o), dok u toploj polovini godine, narocito u junu i julu, u 
Beogradu najvecu cestinu imaju vetrovi iz zapadnog pravca (Unkasevic, 1994). 


Tabela 1. Srednje mesecne destine pravaca vetra i tisina-C (%o) u Beogradu-Opservatorija u 
periodu od 1888. do 1991. godine. 
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Tabela 2. Srednje mesecne brzine vetra (m/s) u Beogradu-Opservatorija u periodu od 1888. do 
1991. godine. 
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Karakteristicno je da se najjaci srednji vetar javlja u martu i to sajacinom 
od 5,9 m/s iz pravca ESE i sajacinom od 4,8 m/s iz SE pravca (tabela 2). Mesec 
mart se odlikuje i sa najvecom cestinom kosavskih perioda u trajanju od 2, 3 i 4 
dana, ali su zabelezeni i periodi u trajanju od 18, 21 i 22 dana (Milosavljevic, 
1950). 


U sinoptickim situacijama kada dolazi do formiranja gradijenta pritiska u 
pravcu od istocne ili jugoistocne Evrope prema zapadnom Sredozemlju, u 
istocnom i severoistocnom delu teritorije Srbije pocinje da duva vetar 
jugoistocnog smera, kosava. To su u vecini slucajeva situacije kada se razvija 
ciklon u zapadnom Sredozemlju (slika 3), ili se nalazi centar anticiklona iznad 
istocnog dela Evrope (slika 4). U prvom slucaju, kosava ima vise juzni smer, 
fenskog je karaktera i zove se topla kosava. U poslednjem slucaju, kosava ima 
vise istocni smer, burskog je karaktera i zove se hladna kosava. 

Poznato je da je fen vetar koji se cesto javlja u Austriji, Svajcarskoj i 
Nemackoj, tj., u oblastima severno od Alpa. To je vetar koji se spusta sa Alpa, 
pri cemu se zagreva i donosi suvo vreme. Posto je kosava u prvom slucaju 
vezana za priblizavanje ciklona sa zapada, ona donosi vazduh sa nizih 
geografskih sirina, tj., ima karakteristike suvog i toplog vetra. Posle ovakve 
kosave obicno dolazi do prodora hladnog vazduha i pojave kise. 

Kada u toku zime dode do formiranja izrazitog anticiklona sa centrom u 
Sibiru, vazduh koji se prebacuje preko planina istocne Srbije ne moze se 
znacajno dinamicki zagrejati i u oblasti Pomoravlja, Podunavlja i istocnog 
Banata duva veoma hladan i slapovit vetar, tipa bure. I u ovom slucaju kosava ne 
donosi padavine osim retkog i suvog snega. 

Kosavu cini vazduh koji se prebacuje preko planina istocne Srbije, a 
narocito je izrazita u oblasti Dunava od Velikog Gradista do Novog Sada i u 
Banatu u okolini Vrsca (Radinovic, 1981). Proucavanjem kosave u Podunavlju 
bavili su se Milosavljevici (1950, 1950a, 1972, 1975, 1976), a kosavom kao 
modifikatorom klime Stanojevic (1959). Ucestanost kosave, odnosno broj dana 
sa kosavom u pojedinim mesecima u hladnoj polovini godine je dosta razlicit od 
godine do godine. Srednji broj dana sa kosavom je u oktobru 8,0, u novembru 
11,3, u decembru 9,1, u januaru 7,3 i u martu 9,1. Prema tome, kosava je 
najcesca u novembru, decembru i martu, a najreda u januaru (Milosavljevic, 
1950a, 1976). 

Trajanje perioda sa kosavom je takode razlicito i zavisi od postojanosti 
sinopticke situacije koja je uslovljava. Najcesce duva 2 do 3 dana, a ponekad i 
znatno duze i do mesec dana. Prosecno traje 5 dana. Prema Milosavljevicu 
(1950, 1976), neprekidno trajanje kosave od 5 dana i duze u periodu od 1920. do 
1945. godine bilo je u oktobru 13 puta, u novembru 21, u decembru 19, u 
januaru 1 1 i u februaru i martu 14 puta. Kosava moze da traje i mesec dana, kada 
pomenuti termicki anticiklon sa centrom iznad Rusije postane blokirajuci. 
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850 hPa Temper alur (Grad C) 


500 hPa Geopotential (gpdm) und Bodendruck (hPa) 


Slika 3. Sinopticke karte Evrope za 14. januar 1972. godine u 01 cas po SEV. 
lokaciju Beograda: a) prizemna karta, b) temperatura na 850 hPa, c) geopotencij 



850 hPa Temperatur (Grad C) 


500 hPa Geopotential (gpdm) und Bodendruck (hPa) 
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Najduzi period sa kosavom u trajanju od 31 dana je zabelezen 1953. godine 
(11. oktobar-10. novembar) i 1972. godine (14. januara-13. februara). Po jacini i 
posledicama po privredu znacajno je pomenuti kosavu od 25. do 29. marta 1959. 
godine, kada su u Vrscu zabelezeni udari vetra od 37,3 m/s, a u Beogradu od 
33, 8 m/s. 

Zbog moguceg uticaja na gradevinske objekte, vise autora se bavilo 
proucavanjem udarnosti vetra bez obzira na pravac. Davenport (1961), 
Borisenko (1977), Wieringa (1973) i Jackson (1978) bavili su se proucavanjem 
udarnosti vetra u prizemnom sloju atmosfere, dok su Davies i Newstein (1968) 
kao i Brook i Spillane (1970) proucavali promenu maksimalnih udara vetra sa 
visinom. Dosadasnja izucavanja udarnosti kosave za teritoriju Srbije odnosila su 
se samo na klimatoloska razmatranja ove pojave (Milosavljevic, 1972; 
Milosavljevic, 1975 i Unkasevic i saradnici, 1998). 

Vukmirovic (1985a) je pokazao da kosava moze biti kako frontalna tako i 
nefrontalna. Proucavajuci prostornu strukturu kosave, Vukmirovic je ustanovio 
da za vreme pojave kosave najvece prostorne promene trpi prizemno polje 
pritiska zbog uticaja orografije pri postojecoj makrocirkulaciji. U istom radu, 
Vukmirovic je ukazao i na zavisnost izmedu vremena pojave dnevnog 
maksimuma intenziteta niske mlazne struje i maksimalnih prizemnih udara 
vetra. Naime, nivo maksimalnog udara kosave se javlja veoma blizu donje 
granice stabilnog sloja, koji obicno pocinje inverzijom. 

Radosavljevic i Vojinovic-Kljaic (1985), koristeci postupak visestruke 
linearne korelacije, su dobili dijagnosticke jednacine za srednju i maksimalnu 
casovnu brzinu vetra koje se mogu koristiti u operativnoj praksi tokom dana sa 
pojavom kosave. 

D. Vukmirovic i saradnici (1997) su na primeru dugotrajne kosave prikazali 
statisticke karakteristike strukture kosave na gradskoj lokaciji. Za vrednosti 
dnevnih brzina vetra, maksimalnih casovnih brzina i faktora udarnosti uradena je 
korelaciona i spektralna analiza da bi se videlo koji opseg ucestanosti navedenih 
velicina ima najvecu energiju. Dobijeno je da najveca energija srednjih brzina 
vetra ima period od oko 122 h, sto odgovara procesima sinoptickih razmera. Za 
faktor udarnosti maksimum je na oko 4h, sto odgovara lokalnim uticajima. 


2 Osnovne karakteristike vetra kosave 


Iznad severoistocnog i istocnog dela Srbije jaka kosava se po pravilu javlja 
kao posledica istovremenog postojanja izrazitog anticiklona sa centrom iznad 
Ukrajine i Besarabije i ciklona u Mediteranu. Ove pojave imaju sezonski 
karakter i najizrazitije su tokom zime. 

Hladan vazduh, koji nadire iz SE pravca nagomilava se i delimicno 
zaustavlja ispred planinskog masiva Karpata da bi kada dostigne kriticnu masu 
poceo da se prebacuje u Vlasku niziju i dalje krece duz korita Dunava, Timoka i 
Crne Reke. Zato se na teritoriji Srbije zahladenje prvo pojavi u Timockoj 
Krajini. Prebacivanje hladnog vazduha preko Homoljskih planina uslovljava 
zatim pojavu kosave u Pomoravlju i Podunavlju. 

Pravac duvanja kosave je uglavnom uslovljen reljefom. Tako u dolini 
Velike Morave i Dunava kosava ima SE pravac dok prema periferiji divergira. U 
ovom slucaju su karakteristicni veliki gradijenti pritiska u podnozju Karpata i 
Homoljskih planina i sve rede izobare (linije istih vrednosti pritiska) prema 
istoku i zapadu. 

Proucavajuci kosavu tokom perioda 1930-1940. godine, Milosavljevic 
(1950) je dosao do zakljucka da prilikom pojave kosave deo vazdusne mase koja 
struji uz Dunav (slika 5, smer a) se odvaja pre Negotina (slika 5, smer b ) i 
prebacuje kroz planinske klance Svrljiskih planina ka Nisu, gde se javlja kao 
fenski padajuci i slapovit vetar velike brzine iz severoistocnog ili istocnog 
pravca (slika 5, smer Z? j ). Preostali deo vazdusne mase koja se kretala uz Timok 
i koja se odvojila od glavne struje i dospela do padina Rtnja, prebacivala se kroz 
planinske klance ka Aleksincu (slika 5, smer Z? 2 ) i ka Paracinu i Cupriji (slika 5, 
smer b 3 ). Za pomenuti period tokom hladne polovine godine (oktobar-mart) 
Milosavljevic je odredio i rezultante pravca vetra za 18 meteoroloskih stanica 
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tokom kosavskih perioda (slika 6). Pravac ovih rezultanti predstavlja pravac 
preovladujuceg vetra na stanici za vreme kosavske situacije, a njihova duzina 
predstavlja cestinu preovladujuceg vetra u %. Poprecne crte na rezultantama 
vetra predstavljaju srednju brzinu kosavskog vetra u m/s, i to duza crta oznacava 
2 m/s, a kraca 1 m/s. To su srednje brzine vetra racunate u periodu 1930-1940. 
godine za najcesce vetrove iz jugoistocnog kvadranta, osim za Nis gde su 
uzimani pravci vetra od severoistocnog do istocnog-jugoistocnog pravca. Kao 
sto se i ocekivalo, najvece brzine kosave u hladnom periodu su registrovane u 
Velikom Gradistu, Beogradu, juznom Banatu i Vrscu. 



Slika 5. Put kosave (Milosavljevic, 1950). 


2 Osnovne karakteristike vetra kosave 
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Slika 6. Rezultante pravca vetra za 18 meteoroloskih stanica tokom kosavskih perioda 

(Milosavljevic, 1950). 

Kao po pravilu, kosava je najcesca u hladnoj polovini godine tj., od oktobra 
do marta, mada se moze javiti i u toploj polovini godine (tabela 3). 
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Kao sto se moze videti iz tabele 3, u toku ispitivanog perioda sa najcescom 
pojavom kosave (1971-1975), broj dana sa kosavom se kretao od 11,5% (leto) 
do 16,1% (zima). Najveca brzina kosave izmerena je tokom jeseni i zime (od 
25,4 do 31,0 m/s), dok je tokom proleca i jeseni maksimalna brzina manja (24,0 
m/s). Pri odabiranju dana sa kosavom Unkasevic i saradnici (1998) su usvojili 
sledeci kriterijum: kao dane u kojima je duvala kosava uzimani su svi oni dani u 
kojima je u sva tri klimatoloska termina osmatranja (07, 14 i 21 cas) duvala 
kosava. Sto se tice intenziteta brzine vetra, u obzir su uzimani samo oni dani u 
kojima je srednja casovna brzina vetra bila jednaka ili veca od 1,0 m/s. 


Tabela 3. Broj dana sa kosavom N (%) i njena maksimalna brzina (m/s) po sezonama u 
Beogradu-Opservatorija u periodu od 1971. do 1975. godine (Unkasevic i saradnici, 1998). 


God. 

Prolede 

N(%) 

v 

max 

Leto 
N(%) l 

v 

r max 

Topla pol. 
god. N(%) 

Jesen 

N(%) 

v 

K max 

Zima 

N(%) 

v 

max 

Hladna pol 
god. N(%) 

1971 

15,7 

24,0 

10,7 ! 

14,8 

26,4 

10,1 

25,4 

10,4 

EK9 

20,5 

1972 

15,7 

19,9 

8,1 


23,9 

11,2 

17,4 

20,5 

mnm 

31,7 

1973 

12,1 

■l&l 

12,4 

1 23,5 | 

24,2 

16,6 

18,5 

13,8 

\ 29,7 1 

30,3 

1974 

19,4 

19,6 

15,4 

mm 

34,8 

19,9 

21,5 

18,8 

EX1 

38,7 

1975 

15,4 


10,7 

16,2 

26,1 

16,6 

23,9 

16,9 

27,6 

33,4 

Sr. vr. 

15,7 


11,5 


27,1 

14,9 


16,1 


30,9 


Treba napomenuti da se najvece brzine vetra u Beogradu javljaju iz 
kosavskog pravca i da je srednja brzina kosave u svim mesecima tokom godine 
skoro dva puta veca od srednje brzine vetra racunate uzimajuci u obzir sve 
pravce duvanja vetra (Stanojevic, 1959; Unkasevic, 2003). 

Analizirajuci kosavu u Beogradu tokom perioda 1920-1945. godine, 
Milosavljavic (1950) je dosao do zakljucka da se kosava u trajanju od pet i vise 
dana najcesce javlja u novembru i decembru, tj. da kosava najduze traje krajem 
jeseni i pocetkom zime. 

Prema srednjim vrednostima broja dana sa kosavom u hladnom periodu 
izgleda da se u novembru javlja duplo vise dana sa kosavom nego u januaru. Iz 
tabele 4 se vidi da je verovatnoca pojave kosave u hladnom periodu (oktobar- 
mart) najveca u trecoj dekadi oktobra i da iznosi 35% (Milosavljevic, 1950a). 
Novembar i decembar imaju najvecu verovatnocu pojave kosave u prvoj, a 
najmanju u trecoj dekadi. Januar, medutim, ima najvecu verovatnocu pojave u 
trecoj dekadi, a najmanju u prvoj. Iz podataka se vidi da je najveca verovatnoca 
pojave kosave u februaru u poslednjih osam dana (31%). U trecoj dekadi marta 
verovatnoca pojave kosave je najveca i iznosi 39%, dok prethodne dve dekade 
imaju skoro upola manju verovatnocu. 


Tabela 4. Verovatnoca pojave kosave u Beogradu za period 1920-1940. god. u procentima (%) 
(Milosavljevic, 1950a). 


Dan/mesec 

X 

XI 

XII 

I 

II 

III 

1 

19 

23 

23 

15 

25 

19 

2 

19 

23 

35 

19 

35 

19 

3 

15 

42 

46 

15 

12 

15 

4 

23 

46 

42 

23 

19 

23 

5 

19 

42 

25 

19 

25 

19 

6 

23 

50 

35 

19 

23 

23 

7 

19 

46 

38 

19 

8 

19 

8 

19 

35 

35 

25 

23 

19 

9 

12 

46 

38 

23 

19 

12 

10 

4 

54 

42 

19 

12 

5 

I dekada 

17 

41 

36 

20 

20 

17 

11 

12 

46 

38 

25 

12 


12 

19 

46 

35 

25 


la 

13 

19 

46 

35 

15 


19 

14 

19 

50 

31 

23 


19 

15 

15 

35 

31 

31 


15 

16 

15 

19 

23 

31 


15 

17 

19 

35 

25 

23 


19 

18 

19 

31 

31 

31 


19 

19 

15 

35 

38 

31 

31 

15 

20 

19 

35 

38 

15 

23 

19 

II dekada 

17 

38 

33 

25 

19 

17 

21 

31 

35 

31 

12 

31 

31 

22 

42 

42 

25 

15 

35 

42 

23 

50 

46 

31 

19 

31 

50 

24 

25 

46 

31 

31 

38 

25 

25 

42 

25 

19 

23 

31 

42 

26 

31 

12 

15 

31 

35 

31 

27 

31 

25 

23 

31 

35 

31 

28 

42 

25 

23 

15 

15 

42 

29 

38 

35 

35 

23 

/ 

38 

30 

25 

31 

35 

31 

/ 

25 

31 

35 

/ 

19 

42 

/ 

35 

III dekada 

35 

32 

26 

27 

31 

39 


Maksimalna verovatnoca racunata po danima (tabela 4) u iznosu od 50% u 
oktobru i martu javlja se 23. dana, sto znaci da ce svake druge godine 23. marta i 
23. oktobra duvati kosava. U novembru je maksimalna verovatnoca desetog, a u 
decembru treceg dana. Maksimum verovatnoce pojave kosave javlja se 
poslednjeg dana januara dok je u februaru osamnaestog dana. 
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U tabeli 5 pokazana su odstupanja pojedinih klimatskih elemenata od 
normalnih vrednosti pri duvanju kosave u periodu od 1920. do 1945. godine 
(Milosavljevic, 1950a). Kao sto se moze videti, najhladnije je pri kosavi tipa 
bure tokom proleca i jeseni, a najtoplije je pri kosavi tipa fena tokom leta, kada 
je osmotren i najveci porast apsolutne vlaznosti i najveci pad relativne vlaznosti. 
Tokom zime za vreme kosave srednje mesecne temperature su vise od 
prosecnih, kao posledica dinamickog zagrevanja vazduha prilikom spustanja, 
dok je relativna vlaznost niza od prosecne. Kosava je padajuci vetar i pri njenoj 
pojavi u istocnim i severoistocnim delovima Srbije vreme je pretezno vedro i 
najveci pad oblacnosti zabelezen je pri kosavi tipa bure tokom leta i jeseni. 
Ukoliko je pojava kosave pracena i izrazitim frontom, oblacnost moze biti 
znacajna kao i izmerena kolicina padavina (slucaj pojave kosave od 16. do 18. 
oktobra 1976. godine, koja je opisana u Glavi 3). 

Interesantno je napomenuti da je leti vetar iz ESE pravca topliji u odnosu na 
vetar iz SE pravca, dok zimi vazi obrnuto (tabela 5). Razlog tome lezi u cinjenici 
sto leti vetar iz ESE pravca na podrucje Srbije donosi tople kontinentalne 
vazdusne mase iz Ukrajinske ravnice, dok vetar iz SE pravca donosi vazdusne 
mase anticiklonskog porekla koje su se formirale iznad Male Azije, Balkana i 
Srednje Azije (oblasti gde se nalaze Sredozemno, Egejsko, Mramorno i Crno 
more). 

Tabela 5. Odstupanje temperature Ar(°C), apsolutne vlaznosti Ae (grrf 3 ), relativne vlaznosti 
A U (%) i oblacnosti An (%) od normalnih vrednosti pri duvanju kosave u Beogradu po 
sezonama (Milosavljevic, 1950a). 


Pravac 

kosave 

Zima 

Prolece 


AT 

Ae 

A U 

An 

AT 

Ae 

A U 

An 

ESE 

0,3 

-0,2 

-4 

-4 

-0,1 

-0,1 

-J 

-1 

SE 

1,9 

0,3 

-6 

-10 

U 

0,4 

-3 

-1 

SSE 

3,9 

0,6 

-9 

-6 

3,3 

0,3 

-9 
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Pravac 

kosave 

Leto 

Jesen 


AT 

Ae 

A U 

An 

AT 


A U 

An 

ESE 

2,6 

0,5 

-3 

-1 

-0,1 

-0,4 

-5 

-14 

SE 

2,1 

0,1 

n 

-13 

1,5 

0,2 

-5 

-9 

SSE 

WSM 



-1 

EH 



-6 


3 Prostorna mezoanaliza kosavskih 
situacija 

U ovoj Glavi monografije analizirana su dva karakteristicna slucaja pojave 
kosave na teritoriji Srbije sa jakim udarima vetra od preko 30,0 m/s i to u 
periodima od 16. do 18. oktobra 1976. godine i od 05. do 06. decembra 1995. 
godine. U prvom slucaju, pojava kosave je bila pracena obilnim padavinama 
(frontalna kosava), dok je za drugi slucaj bilo karakteristicno uglavnom suvo 
vreme (nefrontalna kosava). Kao osnovni materijal za pracenje promena 
osnovnih meteoroloskih elemenata iz casa u cas posluzili su podaci iz obrasca 
“Razvoj vremena” registrovani na glavnim meteoroloskim stanicama lociranim 
u kosavskoj oblasti. U oba slucaja su podaci uzimani sa 18 glavnih 
meteoroloskih stanica i to u prvom slucaju koriscene su stanice: Sombor, Palic, 
Novi Sad, Kikinda, Zrenjanin, Vrsac, Sremska Mitrovica, Beograd, Valjevo, 
Kraljevo, Smederevska Palanka, Veliko Gradiste, Cuprija, Krusevac, Nis, Crni 
Vrh, Negotin i Zajecar, a u drugom slucaju stanice: Palic, Becej, Kikinda, 
Zrenjanin, Vrsac, Banatski Karlovac, Rimski Sancevi, Beograd, Smederevska 
Palanka, Veliko Gradiste, Kragujevac, Cuprija, Krusevac, Crni Vrh, Negotin, 
Nis, Dimitrovgrad i Zajecar. 


3.1 Sinopticka situacija iznad Evrope u periodu 
od 16. do 18. oktobra 1976. godine 

Sinopticka situacija iznad Evrope u periodu od 16. do 18. oktobra 1976. 
godine bila je karakteristicna po sledecem: 
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-16. oktobra 1976. godine u prizemlju (01 cas po SEV) je postojao razvijen 
anticiklon sa centrom u severozapadnom delu Rusije. U Sredozemlju (Tirensko 
more) formirao se mladi ciklon. Teritorija Srbije se nalazila u toplom i vlaznom 
sektoru Sredozemnog ciklona. 

-17. oktobra 1976. godine u prizemlju (01 cas po SEV) Sredozemni ciklon 
je u svom razvoju dostigao nivo od 500 mb na kome postoji zatvorena izohipsa 
(linija koja povezuje mesta iste visine) u oblasti Denovskog zaliva (slika 7a). 
Iznad teritorije Srbije dolazi zbog toga do zgusnjavanja izobara koje su 
orijentisane u pravcu NW-SE, oblacno je i kisovito vreme sa hladnim frontom 
koji se priblizava sa SW (slika 7c). 

-18. oktobra 1976. godine u prizemlju (01 cas po SEV) jaca formiran ciklon 
iznad Velike Britanije, dok ciklon u Sredozemlju slabi uslovljavajuci pad 
gradijenta pritiska (promena pritiska u horizontalnom pravcu) iznad teritorije 
Srbije i smirivanje padavina i kosave. 



Slika 7. Sinopticke katre Evrope: a) AT 500 mb, b) AT 850 mb i c) prizemna karta 
17. oktobra 1976. godine (01 cas SEV). 
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3.2 Polje pritiska iznad kosavske oblasti u 
periodu od 16. do 18. oktobra 1976. godine 

Detaljnu analizu kosavske situacije iznad teritorije Srbije u periodu od 16. 
do 18. oktobra 1976. godine uradilaje Petrovic (1977). Posmatrajuci tendenciju 
pritiska (mb) u pet termina u trajanju od po 12 casova (tabela 6) mozemo da 
zakljucimo sledece: 

-u prvom terminu postojao je pad pritiska na zapadnom i porast pritiska na 
istocnom delu kosavske oblasti; 

-u drugom terminu pritisak je rastao u citavoj kosavskoj oblasti; 

-u trecem terminu postojao je najveci porast pritiska u citavoj kosavskoj 
oblasti; 

-u cetvrtom terminu u istocnom delu kosavske oblasti pritisak je poceo da 
opada, sto je ukazivalo na pocetak slabljenja posmatranog mezosistema; 

-u petom terminu uocen je opsti porast pritiska sa najvecim intenzitetom u 
zapadnom delu kosavske oblasti. 

3.3 Vetar u kosavskoj oblasti u periodu od 16. 
do 18. oktobra 1976. godine 

Posmatrajuci pravac vetra u kosavskoj oblasti u periodu od 16. do 18. 
oktobra 1976. godine (tabela 7) mozemo zakljuciti sledece: 

-Beograd je jedino mesto sa konstantnim pravcem vetra od 16. do 18. 
oktobra, pa se moze smatrati centrom kosavske oblasti; 

-mala odstupanja od srednjeg dnevnog pravca vetra imala su mesta koja 
blize ili dalje okruzuju Beograd (Vrsac, Veliko Gradiste, Zrenjanin, Novi Sad, 
Smederevska Palanka i Kraljevo); 

-najveca odstupanja od srednjeg dnevnog pravca vetra imala su mesta koja 
se nalaze na periferiji kosavske oblasti (Palic, Kikinda, Valjevo, Krusevac, 
Cuprija, Nis i Negotin). 

Najjaci udari kosave iznad teritorije Srbije javljali su se u oblasti Vrsca i 
Velikog Gradista odakle su slabili ka oblastima Podunavlja i Pomoravlja. 
Prostorno maksimalan domet izotaha (linije koje povezuju mesta sa istom 
brzinom vetra) dostignut je u periodu maksimalnog intenziteta vetra (u Vrscu je 


19 


3 Prostorna mezoanaliza kosavskih situacija 

zabelezen maksimum od 35 m/s) izmedu 07 i 19 casova 17. oktobra 1976. 
godine. Sa jacanjem intenziteta vetar je postajao stalniji po pravcu. 
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3.4 Temperatura vazduha u kosavskoj oblasti u 
periodu od 16. do 18. oktobra 1976. godine 

U tabeli 8 prikazane su trocasovne promene temperature vazduha izmerene 
na 18 glavnih meteoroloskih stanica u kosavskom podrucju u periodu od 16. do 
1 8. oktobra 1976. godine. Iz tabele 8 se moze videti da su se tokom posmatranog 
perioda javila dva talasa naglog otopljenja. Jedan je bio 16. oktobra 1976. godine 
u intervalu od 07 do 13 casova i prethodio je naglom padu temperature. Drugi je 
bio 18. oktobra 1976. godine u intervalu od 01 do 13 casova i javio se kao 
posledica prestanka jake kosave. 

Iz gore pomenute tabele moze se uociti da je 16. oktobra 1976. godine u 
intervalu od 13 do 19 casova u vecini mesta kosavske oblasti doslo do naglog 
pada temperature vazduha. Zapazeno je da je trocasovni pad u oblasti Vrsca, 
Velikog Gradista, Cuprije i Negotina bio manji po intenzitetu od pada 
temperature u preostalom delu posmatrane kosavske oblasti. Najveci pad 
temperature je zabelezen 16. oktobra u intervalu od 13 do 16 casova u 
centralnom delu kosavske oblasti, a kasnije u intervalu od 1 6 do 19 casova po 
granici kosavske oblasti. Moze se zakljuciti da je ukupan trodnevni pad 
temperature, kao posledica pojave kosave, bio najveci u oblasti zapadno od 
Velike Morave i juzno od Dunava. Najmanja razlika izmedu maksimalne i 
minimalne temperature vazduha (10°C) uocena je u okolini Vrsca. 

Takode je uoceno da se nagli pad temperature (16. oktobra u periodu od 13 
do 16 casova) u centru kosavske oblasti javljao istovremeno sa nailaskom 
izotahe od 10,0 m/s. Na periferiji kosavske oblasti najveci trocasovni pad 
temperature (16. oktobra u periodu od 16 do 19 casova) je takode bio pracen 
izotahom od 10,0 m/s. Sa prestankom kosave doslo je do porasta temperature u 
citavoj kosavskoj oblasti. 
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Na osnovu radiosondaznih podataka o atmosferi iznad Beograda (u 01 i 13 
casova) u periodu od 16. do 18. oktobra 1976. godine Petrovic (1977) je ukazala 
na sledece: 

-16. oktobra 1976. godine u 01 cas postojala je temperaturska inverzija 
(porast temperature sa visinom) do visine od 820 m, dok se u 13 casova inverzija 
podigla i zauzela sloj visine od 620 do 930 m. 

-17. oktobra 1976. godine u 01 cas (slika 8) vrlo izrazena inverzija 
(5°C/100 m) prostirala se u sloju visine od 600 do 900 m usled prodora hladnog 
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vazduha u plitkom prizemnom sloju i jakog smicanja vetra na nivou od oko 910 
mb dok se u 13 casova inverzija podigla i nalazila se u sloju od 760 do 1370 m, 
sto su pokazali i Vojnovic Kljaic i Popovic (1978). 



Slika 8. Radiosondazno merenje u Beogradu u 01 cas, 17. oktobra 1976. god. (Vojnovic Kljaic i 

Popovic, 1978) 


-18. oktobra 1976. godine temperaturska inverzija intenziteta 3°C/100 m 
prostirala se u sloju visine od 550 do 910 m (SE vetar je registrovan do visine od 
3000 m), dok su se u 13 casova javile dve inverzije: jedna izmedu visina od 550 
do 1 150 m, a druga izmedu visina od 1340 do 1470 m (SE vetar je registrovan 
do visine od 1100 m). 

3.5 Padavine u kosavskoj oblasti u periodu od 
16. do 18. oktobra 1976. godine 

U periodu od 19 casova 16. oktobra 1976. godine do 01 casa 17. oktobra 
1976. godine zabelezene su padavine u citavoj kosavskoj oblasti. 

Sa slike 9 se moze videti da su izohijete (linije koje povezuju iste vrednosti 
padavina) orijentisane u pravcu NE-SW, tj. normalno na pravac osnovnog 
strujanja vazduha. Od SE ka NW kolicina padavina je opadala sve do Kraljeva, 
Cuprije, Vrsca i Velikog Gradista, kao posledica spustanja vazduha. Kolicine 
padavina su zatim rasle prema Zrenjaninu, Novom Sadu i Sremskoj Mitrovici, 
pa se opet smanjivale prema Somboru i Palicu. Slican raspored padavina sa 
naizmenicnim porastima i padovima osmotrio je i Flohn (Radinovic, 1979) u 
vezi sa intertropskim frontom (zamisljena linija koja deli dve vazdusne mase) 
iznad Atlantika, objasnjavajuci to dinamikom polja vetra. 

U periodu od 01 do 07 casova 17. oktobra 1976. godine padavine su bile 
skoncentrisane na oblast Homolja opadajuci ka zapadu, dok je u Pomoravlju i 
Podunavlju bilo suvo vreme sto je ukazivalo na pocetak kulminacije kosave. U 
periodu od 07 do 13 casova 17. oktobra 1976. godine minimum padavina je 
registrovan u oblasti najjaceg vetra (Vrsac, Veliko Gradiste i Beograd) i iznosio 
je 1 mm, dok je iduci prema severu kolicina padavina rasla. 

Sa jacanjem kosave se kolicina padavina smanjivala. Tako je u periodu od 
13 do 19 casova 17. oktobra 1976. godine najmanje padavina bilo u Vrscu (0,2 
mm) i u oblasti Zajecara, Crnog Vrha i Nisa (do 1,0 mm), dok je u Beogradu 
izmereno nesto vise tj., 3,6 mm. 
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Slika 9. Kolicine padavina (mm) u kosavskoj oblasti u periodu od 16. oktobra (19 casova) do 

17. oktobra (01 cas) 1976. godine (Petrovic, 1977). slika 10 - Kolicine padavina (mm) u kosavskoj oblasti izmerene u periodu od 16. do 18. oktobra 

1976. godine (Petrovic, 1977). 
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Na slici 10, na kojoj su prikazane ukupne trodnevne padavine, je evidentno 
da je najmanje padavina (od 5,0 do 10,0 mm) izmereno u oblasti koja obuhvata 
Vrsac, Beograd, Sremsku Mitrovicu, Valjevo, Cupriju i Veliko Gradiste tj., 
oblast sa najjacom kosavom (slika 1 1). 



Slika 11. Izohrone (mesta sa istim casovnim vrednostima brzine) izotahe (mesta sa istim 
vrednostima brzine) 20 in/s (Petrovic, 1977). 


Prema severu i istoku kolicina padavina se povecavala i dostigla je 
maksimalnu vrednost u Negotinu (53,8 mm). Ovakav raspored padavina je 
posledica postojece orografije, dok se porast padavina na severu moze objasniti 
smanjenjem fenskog efekta vetra kosave u postojecoj ciklonalnoj situaciji. 

Tokom analize posmatrane kosavske situacije od 16. do 18. oktobra 1976. 
godine zapazeno je da su oblasti sa konstantnim pravcem vetra istovremeno i 
oblasti sa najjacim vetrom, sto ukazuje da ce u toj oblasti biti izmerene najmanje 
padavine. Medutim, 17. oktobra 1976. godine izmerene su najvece padavine 
tokom trodnevnog perioda pojave kosave, kada je bila kulminacija sistema. 
Posto su u posmatranoj kosavskoj situaciji analizirani maksimalni casovni udari 
vetra koji su kratkog trajanja, to su se padavine 17. oktobra 1976. godine 
verovatno izlucile u periodima slabijeg vetra izmedu maksimalnih udara. 

3.6 Sinopticka situacija iznad Evrope u periodu 
od 05. do 06. decembra 1995. godine 

Sinopticka situacija iznad Evrope u periodu od 05. do 06. decembra 1995. 
godine bila je karakteristicna po sledecem: 

-05. decembra 1995. godine u prizemlju je postojao veoma razvijen 
anticiklon sa centrom iznad Skandinavije (1045 mb u centru). Jugoistocno od 
nase zemlje postojao je slab ciklon sa centrom u Egejskom moru, dok se u 
zapadnom Sredozemlju formirao mladi ciklon (slika 12). Teritorija Srbije se 
nalazila u relativno toplom i vlaznom vazduhu, kada je preovladivalo oblacno 
vreme praceno slabim snegom. Tokom sledecih 12 casova anticiklon se iz 
oblasti Skandinavije prosirio i pomerio na SE obuhvatajuci tako celu Evropu 
osim Sredozemlja i veci deo bivseg Sovjetskog Saveza. Ciklon u Egejskom 
moru se ugasio, dok se ciklon iz zapadnog premestio u oblast centralnog 
Sredozemlja dobijajuci na intenzitetu (zatvorene izohipse su ocigledne na 500 
mb nivou). 

-06. decembra 1995. godine u prizemlju je postojao razvijen ciklon sa 
centrom u Tirenskom moru (slika 13). Iznad teritorije nase zemlje izobare su 
bile orijentisane u pravcu NW-SE, a preovladivalo je suvo i oblacno vreme. 
Tokom sledecih 12 casova ciklon u centralnom Sredozemlju je oslabio, sto je 
imalo za posledicu pad gradijenta pritiska iznad teritorije Srbije i slabljenje 
kosave. 
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Slika 12. Prizemna sinopticka situacija iznad Evrope 05. decembra 1995. godine 

(00 casova SEV). 


00 


Slika 13. Prizemna sinopticka situacija iznad Evrope 06. decembra 1995. godine 

(00 casova SEV). 
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3.7 Vetar u kosavskoj oblasti u periodu od 05. 
do 06. decembra 1995. godine 

Posmatrajuci pravac vetra tokom dva dana pojave kosave (tabela 9) 
mozemo zakljuciti da je: 

-Beograd jedino mesto sa konstantnim pravcem vetra u periodu od 05. do 
06. decembra 1995. godine, pa se kao i u prethodnom slucaju pojave kosave, 
moze smatrati centrom kosavske oblasti; 

-najmanja odstupanja od srednjeg dnevnog pravca vetra imala su mesta koja 
se nalaze najblize Beogradu (Rimski Sancevi, Banatski Karlovac, Zrenjanin, 
Becej, Kikinda,Vrsac,Veliko Gradiste i Smederevska Palanka); 

-najveca odstupanja od srednjeg dnevnog pravca vetra imala su mesta koja 
se nalaze na periferiji kosavske oblasti (Palic, Kragujevac, Cuprija, Krusevac, 
Negotin, Zajecar, Nis i Dimitrovgrad). 

Interesantno je napomenuti da se drugog dana kosava prosirila dopiruci na 
sever do Palica, a na jug do Cuprije. Istovremeno sa zahvatanjem vece povrsine, 
kosava je dobila na intenzitetu tako da se srednja brzina vetra iz kosavskog 
pravca 06. decembra 1995. godine kretala od 5,5 m/s u Beceju do 11,5 m/s 
koliko je izmereno u Vrscu i Velikom Gradistu. 

Najjaci udari kosave iznad teritorije Srbije zabelezeni su 06. decembra 
1995. godine u trouglu izmedu Vrsca, Beograda i Velikog Gradista odakle su 
slabili ka oblastima Pomoravlja i Podunavlja. Drugog dana duvanja kosave sa 
jacanjem intenziteta, vetar je postajao stalniji po pravcu (tabela 9). Kao i u 
slucaju frontalne kosave, najjaci udari su iznosili 28 m/s i izmereni su 06. 
decembra izmedu 16 i 17 casova u Vrscu. 

3.8 Temperatura vazduha u kosavskoj oblasti u 
periodu od 05. do 06. decembra 1995. 
godine 

U tabeli 10 prikazane su sest-casovne promene temperature vazduha 
izmerene na 18 meteoroloskih stanica u kosavskom podrucju u periodu od 05. 
do 06. decembra 1995. godine. 
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Kao sto se iz tabele 10 moze videti, na pocetku duvanja kosave 05. 
decembra 1995. godine u periodu od 01 do 07 casova u svim mestima u 
kosavskoj oblasti temperatura je opala. Zapaza se da je sest-casovni pad 
temperature u oblasti Kikinde, Banatskog Karlovca, Vrsca, Velikog Gradista, 
Smederevske Palanke, Krusevca, Crnog Vrha, Negotina i Zajecara bio veci po 
intenzitetu od pada temperature u preostalom delu kosavske oblasti. 

Sa jacanjem intenziteta vetra, doslo je do postepenog porasta temperature i 
to prvo 05. decembra 1995. godine u periodu od 13 do 19 casova u Beceju, 
Zrenjaninu, Rimskim Sancevima i Beogradu, azatim 06. decembra 1995. godine 
u Rimskim Sancevima, Banatskom Karlovcu, Vrscu, Nisu i Dimitrovgradu u 
periodu od 01 do 07 casova. Sestog decembra 1995. godine u periodu od 07 do 
13 casova u svim mestima u kosavskoj oblasti doslo je do porasta temperature 
vazduha. Interesantno je da je sest-casovni porast temperature vazduha bio manji 
po intenzitetu u centru nego na obodu kosavske oblasti. Navedena cinjenica se 
moze objasniti porastom srednjih brzina vetra u centru kosavske oblasti, koje su 
se 06. decembra 1995. godine u pomenutom vremenskom periodu kretale u 
intervalu od 7 do 13 m/s (slika 14). Porast temperature se zadrzao istog dana i u 
vremenskom periodu od 13 do 19 casova i to u Palicu, Kikindi, Vrscu i Velikom 
Gradistu, dok je u preostalom, vecem delu kosavske oblasti temperatura vazduha 
bila u padu. 

Za razliku od prethodnog slucaja, kada smo posmatrali frontalnu i hladnu 
kosavu (tipa bure), ovaj slucaj pojave kosave bio je karakteristican po odsustvu 
padavina i porastu temperature vazduha u kosavskoj oblasti. To znaci da smo 
analizirali nefrontalnu kosavu (tipa fena) koja je bila uslovljena adijabatskim 
spustanjem vazduha preko Karpata ka podrucju istocne Srbije. 


3.9 Relativna vlaznost vazduha u kosavskoj 
oblasti u periodu od 05. do 06. decembra 
1995. godine 

U tabeli 1 1 prikazane su sest-casovne promene relativne vlaznosti izmerene 
na 18 meteoroloskih stanica u kosavskom podrucju u periodu od 05. do 06. 
decembra 1995. godine. 

Iz pomenute tabele se moze videti da je prvog dana pojave kosave (05. 
decembar 1995. godine) zabelezen pad relativne vlaznosti skoro u celoj 
kosavskoj oblasti. U vremenskom periodu od 07 do 13 casova najveci pad 
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relativne vlaznosti (od 13 do 21%) zabelezen je u mestima sa izrazitim porastom 
temperature vazduha (Vrsac, Smederevska Palanka, Kragujevac, Cuprija i 
Krusevac) i umerenom kosavom (od 3 do 8 m/s). Istog dana, u vremenskom 
intervalu od 13 do 19 casova pad relativne vlaznosti je bio izrazen u Beceju, 
Zrenjaninu, Banatskom Karlovcu, Beogradu i Velikom Gradistu. 
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Slika 14. Srednje casovne brzine (tanka linija) i pravac vetra (deblja linija) u kosavskoj oblasti 

u 13 casova 06. decembra 1995. god. 
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Interesantno je napomenuti da se sa jacanjem kosave 06. decembra 1995. 
godine pad relativne vlaznosti smanjivao u celoj kosavskoj oblasti i u periodu od 
13 do 19 casova samo je na teritoriji Vojvodine registrovan suvlji vazduh, dok je 
u ostalim mestima kosavske oblasti zabelezen porast relativne vlaznosti (od 2 do 
7%). 
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4 Statisticka analiza maksimalnih 
brzina i faktora udarnosti vetra 
kosave u Beogradu 


4.1 Korisceni podaci 

Za analizu maksimalnih brzina i faktora udarnosti vetra kosave korisceni su 
podaci vetra izmereni na dve meteoroloske stanice: Beograd-Opservatorija i 
Beograd-Zeleno Brdo u periodu najcesce pojave vetra kosave tj od 1971. do 
1975. godine. Kao sto je poznato, meteoroloska stanica Beograd-Opservatorija 
locirana je u uzem centru grada (sa visinom anemometra od 24 m, oznacena sa 1 
na slici 15), a stanica Beograd-Zeleno Brdo nalazila se oko 6 km jugoistocno od 
centra grada (77=243 m sa visinom anemometra od 1 8 m, oznacena sa 2 na slici 
15). Treba napomenuti da je meteoroloska stanica Beograd-Opservatorija tipicna 
gradska stanica, dok se stanica Beograd-Zeleno-Brdo u pomenutom petogo- 
disnjem periodu ispitivanja nalazila na samom obodu grada pa je imala karakter 
prigradske stanice. 

Meteoroloska stanica Beograd-Opservatorija je sa zapada, severa i istoka 
zaklonjena visokim zgradama, dok je prema jugoistoku skoro potpuno otvorena. 
Pod uticajem okolnih zgrada vetar iz kosavskog pravca dobija uzlaznu 
komponentu, cime gubi na jacini. To znaci da su srednje brzine vetra kosave 
izmerene na stanici Beograd-Opservatorija manje nego na stanici Beograd- 
Zeleno Brdo, koja se nalazila skoro u neporemecenoj atmosferi i na nadmorskoj 
visini vecoj za 1 1 1 m u odnosu na stanicu Beograd-Opservatorija. 
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Slika 15. Lokacije meteoroloskih stanica u Beogradu: Meteoroloska Opservatorija na Vracaru 
0 i Aeroloska Opservatorija na Zelenom Brdu 0. 


4.2 Statisticka analiza maksimalnih brzina vetra 
kosave 

Statisticka analiza maksimalnih dnevnih brzina vetra u danima sa pojavom 
kosave u Beogradu u periodu od 1971. do 1975. godine radena je primenom 
sledecih funkcija gustine raspodele (FGR): 

(a) f\x) = ^=expf - ^ nx \ dvoparametarske lognormalne raspodele 


(a) /(*) = 7== ex P — y 

xc-j2n ^ 2 a 

sa parametrima pi a ; 


(b) /(*) = — exp - , troparametarske Vejbulove (Weibull) 

b\ b ) \\ b )) 

raspodele sa parametrima a,b i c ; 

a « ( x y~ ] ( ax' 

(c) f(x) = — exp , dvoparametarske Gama (Gamma) 

raspodele sa parametrima a ip; 

(d) /(*) = a exp {- a(x - p) - exp[- a(x - P) ] } , dvoparametarske Gumbelove 
(Gumbel) raspodele sa parametrima a i P . 

Dobijeni rezultati za meteoroloske stanice Beograd-Zeleno Brdo i Beograd- 
Opservatorija, tj., empirijski i teorijski % 2 n testovi za koriscene funkcije gustine 


dvoparametarske Gama (Gamma) 
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raspodele (sa a =0,05 stepeni slobode) pokazani su u tabeli 12 (Malisic i 
Unkasevic, 1997). 


2 2 

Tabela 12. Empirijski % n i teorijski X w>a =o ,05 statisticki testovi za odabrane funkcije gustine 
raspodele u danima sa kosavom u Beogradu u periodu od 1971 . do 1975. godine. 


Teorijske funkcije 

Beograd-Zeleno Brdo 

Beograd-C 

Ipservatorija 

gustine raspodele 

Empirijska 

Teorijska 

Empirijska 

Teorijska 

fix) 

xl statistika 

Xw;o.os statistika 

xl statistika 

x'.o.os statistika 

Lognormalna 

31,754 

38,885 

62,354 

40,113 

Vejbulova 

163,796 

37,652 

25,405 

38,885 

Gama 

45,924 

38,885 

35,080 

40,113 

Gumbelova 

32,603 

38,885 

57,075 

40,113 


Kao sto se moze videti iz tabele 12 i sa slike 16 dvoparametarska Iog- 
normalna FGR pokazuje najbolje slaganje sa empirijskim podacima izmerenim 
na stanici Beograd-Zeleno Brdo. Medutim, ustanovljeno je da za stanicu 
Beograd-Opservatorija najbolje slaganje sa izmerenim podacima o trenutnoj 
maksimalnoj brzini vetra pokazuje troparametarska Vejbulova FGR (slika 17, 
Unkasevic i saradnici, 2002). Prema tabeli 12 uslov xl -xf, a = oos zadovoljavaju 

takode i Gumbelova FGR za podatke izmerene na stanici Beograd-Opservatorija 
kao i Gama FGR za podatke izmerene na stanici Beograd-Zeleno Brdo. Uoceno 
neslaganje u opisivanju trenutnih maksimalnih brzina vetra iz kosavskog pravca 
posledica je urbanog uticaja na strujanje vazduha na mestu Beograd- 
Opservatorija. Naime, meteoroloska stanica Beograd-Zeleno Brdo nalazila se u 
proucavanom periodu (1971-1975) na samoj granici grada i na znatno vecoj 
nadmorskoj visini u odnosu na stanicu Beograd-Opservatorija. To znaci da se za 
opisivanje trenutnih maksimalnih brzina vetra iz kosavskog pravca u urbanim 
sredinama moze primeniti troparametarska Vejbulova FGR, dok se u slucaju 
skoro neporemecene atmosfere moze uspesno primeniti dvoparametarska 
lognormalna FGR, sto je ustanovljeno i kod drugih autora (Takle and Brown, 
1978; Labraga, 1994; Nof and Paldor, 1994). 
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Slika 16. Poredenje dvoparametarske lognormalne sa: a) Gama, b) Gumbelovom i 
c) troparametarskom Vejbulovom FGR i empirijskim histogramima relativne ucestanosti 
trenutnih maksimalnih brzina vetra izmerenih na stanici Beograd-Zeleno Brdo u periodu od 

1971. do 1975. godine. 


Slika 17. Poredenje troparametarske Vejbulove sa: a) Gama, b) Gumbelovom i 
c) dvoparametarskom lognormalnom FGR i empirijskim histogramima relativne ucestanosti 
trenutnih maksimalnih brzina vetra izmerenih na stanici Beograd-Opservatorija u periodu od 

1971. do 1975. godine. 








40 


Unkasevic, Tosic: Kosava 


4 Statisticka analiza kosave 


41 


Najvece relativne ucestanosti trenutnih maksimalnih brzina kosave za 
stanicu Beograd-Zeleno Brdo osmotrene su u intervalu brzina od 8,5 do 9,5 m/s 
(slika 16), a za stanicu Beograd-Opservatorija u intervalu brzina od 7,5 do 8,5 
m/s (slika 17). Medutim, srednje casovne brzine vetra pokazuju najvece 
relativne ucestanosti u intervalu brzina od 3,0 do 4,0 m/s (slika 18). 

Takode je ustanovljeno da se i srednje casovne brzine vetra, u casu 
osmotrenih trenutnih maksimalnih brzina vetra iz kosavskog pravca na stanici 
Beograd-Opservatorija, mogu dobro opisati pomocu troparametarske Vejbulove 

raspodele (slika 18). U ovom slucaju je takode zadovoljen uslov y jn — XAi,a=o,05 

tj., 12,740 < 16,919. 

Empirijske i teorijske statistike vetra kosave po sezonama u Beogradu- 
Opservatorija (1971-1975) predstavljene su u tabeli 13. 

Primenjujuci vec navedene cetiri funkcije gustine raspodele dobijeno je da 
se dvoparametarska lognormalna FGR najbolje prilagodava trenutnim 
maksimalnim brzinama vetra u svim sezonama izmerenim na stanici Beograd- 
Zeleno Brdo (slika 19). Medutim, na stanici Beograd-Opservatorija utvrdeno je 
da se, u toploj polovini godine, empirijskim podacima najbolje prilagodava 
Gama, a u hladnoj polovini godine troparametarska Vejbulova FGR (slika 20). 



Slika 18. Troparametarska Vejbulova FGR i histogram relativne ucestanosti srednjih casovnih 
brzina vetra izmerenih na stanici Beograd-Opservatorija u periodu od 1971. do 1975. godine. 






Slika 19. Empirijska i dvoparametarska lognormalna FGR trenutnih maksimalnih brzina vetra 
izmerenih na stanici Beograd-Zeleno Brdo po sezonama u periodu od 1971. do 1975. godine. 







42 


Unkasevic, Tosic: Kosava 




Slika 20. Sezonska empirijska i a), b) Gamma i c), d) Vejbulova FGR trenutnih maksimalnih 
brzina vetra izmerenih na stanici Beograd-Opservatorija u periodu od 1971. do 1975. godine. 
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Za ocekivati je da se u relativno neporemecenoj atmosferi trenutna 
maksimalna brzina kosave tokom cele godine najbolje moze opisati samo 
jednom teorijskom funkcijom gustine raspodele, kao i da u gradskoj atmosferi to 
nije slucaj, sto su izmereni podaci i potvrdili. 


Tabela 13. Sezonski empirijski xl i teorijski a =oos testovi za odabrane funkcije gustine 
raspodele u danima sa kpsavom u Beogradu u periodu od 1971 . do 1975. godine. 



Beograd-Zeleno Brdo 

Beograd-( 

Dpservatorija 

Sezona 

Empirijska 
x] statistika 

Teorijska xl.o.os 
statistika 

Empirijska 
x] statistika 

Teorijska x, 2 l;0 05 
statistika 

Prolece 

Lognormalna 

20,447 

Lognormalna 

31,410 

Gama 

14,303 

Gama 

31,410 

Leto 

Lognormalna 

22,946 

Lognormalna 

30,114 

Gama 

21,511 

Gama 

31,410 

Jesen 

Lognormalna 

11,375 

Lognormalna 

32,671 

Vejbulova 

19,892 

Vejbulova 

30,114 

Zima 

Lognormalna 

36,629 

Lognormalna 

43,773 

Vejbulova 

43,419 

Vejbulova 

38,885 


4.3 Analiza udarnosti vetra kosave 

Poznato je da kosava duva na udare, tj., da je mahovite prirode, odnosno 
turbulentne strukture. To se ogleda u naglim i kratkotrajnim oscilacijama brzine 
i pravca oko neke srednje vrednosti, sto se vidi iz zapisa instrumenta, 
anemografa (slika 1). 

Mikroturbulentna struktura umerenih i jakih vetrova, kao sto je kosava, 
predstavljena je fluktuacijama brzine i pravca vetra. Fluktuacija brzine ili 
udarnosti vetra nastaje kao posledica turbulencije u sloju trenja, koja je cesto 
modifikovana efektima orografije kao sto je slucaj kod pojave kosave (Cupic sa 
saradnicima, 1985). Turbulentna struktura se formira usled superpozicije 
turbulentnih fluktuacija komponenata brzine vetra u, v i w najrazlicitijih 
dimenzija na osnovno strujanje. Medutim, u atmosferi mogu da se jave i krupniji 
vrtlozi duzeg trajanja, tzv., koherentna vrtlozna strujanja (Vukmirovic i Merkle, 
1991). 
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Stepen mikrofluktuacije izrazava se indeksom udarnosti definisanim 
pomocu: 

G = 

V 

gde je v max maksimalna brzina vetra, a v srednja brzina vetra u casu kada je 
osmotrena maksimalna brzina vetra. Na osnovu merenja (Borisenko, 1977) 
ustanovljena je linearna zavisnost izmedu maksimalnih i srednjih casovnih 
brzina vetra u obliku: 

v max =OV+b, 

gde su a i b konstante. Zamenjujuci ovaj izraz u jednacinu za indeks udarnosti 
dobijamo da je: 

G = a + k, 

v 

sto takode ukazuje na linearnu zavisnost. 

Informacija o indeksu udarnosti je neophodan podatak za izracunavanje 
dinamickog dela pritiska vetra na objekte u gradu, a posebno dalekovode. Stoga 
su za utvrdivanje zavisnosti izmedu maksimalne i srednje casovne brzine vetra 
iz kosavskog pravca korisceni podaci merenja na stanici Beograd-Opservatorija 
u istom petogodisnjem periodu (1971-1975) kao i u slucaju ispitivanja ponasanja 
trenutnih maksimalnih brzina vetra. U navedenom periodu uocen je i rekordan 
period sa neprekidnom pojavom kosave u trajanju od 31 dana (14. januar-13. 
februar 1972. godine), pa je ova situacija posebno analizirana i dobijeni rezultati 
su uporedeni sa analognim rezultatima iz petogodisnjeg perioda. 

Empirijska veza izmedu maksimalnih dnevnih i srednjih casovnih brzina 
vetra iz kosavskog pravca u petogodisnjem periodu pokazana je na slici 21a, a 
analogna zavisnost izmedu maksimalnih casovnih i srednjih casovnih brzina 
vetra kosave u trajanju od 31 dana pokazana je na slici 21b. 

Zajednicko za obe situacije je postojanje znacajne korelacije izmedu 
ispitivanih velicina sa koeficijentom korelacije vecim od 0,90, kao i skoro isti 
nagib linije regresije prema apscisnoj osi (slika 21). Navedeni rezultat ukazuje 
mozda na zakljucak da bilo koji slucaj kosave dovoljno dugog trajanja sadrzi u 
sebi sve znacajne karakteristike kosave dobijene ispitivanjem na osnovu 
visegodisnjih osmatranja. 

Merenja brzine vetra na razlicitim mestima u svetu pokazala su da su 
maksimalne brzine vetra na nivou merenja, bez obzira na pravac, proporcionalne 


srednjoj brzini vetra u sloju atmosfere od nekoliko stotina metara iznad nivoa 
merenja (Frankenberger, 1968; Borisenko, 1977). Tako se i maksimalne brzine 
kosave na nivou anemografa javljaju istovremeno sa formiranjem jezgra mlazne 
struje sa osom na visini izmedu 250 i 700 m iznad tla (Vukmirovic, 1985a,b). 




Slika 21 . Zavisnost izmedu: a) maksimalnih dnevnih ( v max ) i srednjih casovnih ( v ) brzina 

kosave u periodu od 1971. do 1975. godine i b) maksimalnih casovnih ( v max ) i srednjih 

casovnih ( v ) brzina kosave u periodu od 14. januara do 13. februara 1972. godine izmerenih na 

stanici Beograd-Opservatorija. 
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Maksimalni udari kosave u ispitivanom petogodisnjem periodu javljali su se 
u razlicito doba dana i noci, sto ukazuje na cirkulacione uslove njenog nastanka 
(Vukmirovic, 1985a,b). Udari vetra su bili jaci kod kosave koju je pratio 
prolazak fronta nego kod nefrontalne kosave, kao sto je vec istaknuto u Glavi I 
ove Monografije. 

Takode je ispitivana i relacija izmedu maksimalnih dnevnih i srednjih 
casovnih brzina kosave izmerenih na stanici Beograd-Opservatorija tokom sve 
cetiri sezone u godini, u istom petogodisnjem periodu, sto je pokazano na slici 
22 (Unkasevic i saradnici, 1998). 




izmerenih na stanici Beograd-Opservatorija u razlicitim sezonama u periodu 
od 1971. do 1975. godine. 

Interesantno je da se tokom leta i zime (slika 22b, d) javljaju nagibi 
regresije skoro jednaki nagibu regresije dobijenom za celokupni ispitivani 
petogodisnji period (slika 21a). Razliciti nagibi regresije tokom proleca i jeseni 
(slika 22a, c) verovatno su uslovljeni brzim smenama razlicitih sinoptickih 


situacija iznad Srbije u prelaznim sezonama. Koeficijenti korelacije izmedu 
maksimalnih dnevnih i srednjih casovnih brzina kosave u razlicitim sezonama 
krecu se u intervalu od 0,8927 (leto) do 0,9579 (zima). 

Zavisnost koeficijenta udarnosti G od srednje casovne brzine vetra kosave 
(v) izmerene na stanici Beograd-Opservatorija u ispitivanom petogodisnjem 
periodu pokazana je na slici 23a, a analogna zavisnost za 31-dnevnu kosavu 
pokazana je na slici 23b (Unkasevic i Tosic, 1996). Moze se uociti znacajna 
slicnost oba grafika. Naime, za srednje casovne brzine vetra vece od 12 m/s, 
koeficijent udarnosti G ne zavisi od srednje casovne brzine kosave i ima 
vrednost 2,15, odnosno 2,11, respektivno. 




Slika 23. Zavisnost koeficijenta udarnosti (G ) od srednje casovne brzine vetra kosave ( v ) za 
stanicu Beograd-Opservatorija u periodima: a) 1971-1975. godina, b) 14. januar-13. februar 

1972. godine. 
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Funkcionalna veza izmedu koeficijenta udarnosti i srednje casovne brzine 
vetra kosave izmerene na stanici Beograd-Opservatorija u razlicitim sezonama 
pokazana je na slici 24. Uoceno je da u toploj polovini godine, koeficijenti 
udarnosti imaju vece vrednosti pri malim srednjim casovnim brzinama vetra 
nego u hladnoj polovini godine. Takode, u prolece^ leto, zavisnost koeficijenta 
udarnosti (G ) od srednjih casovnih brzina brzina ( v ) prestaje pri nesto manjim 
vrednostima v nego u jesen i zimu (slika 24). 




Slika 24. Zavisnost koeficijenta udarnosti ( G ) od srednje casovne brzine vetra kosave ( v ) 
izmerene na stanici Beograd-Opservatorija u razlicitim sezonama u periodu od 1971. do 1975. 

godine. 


5 Vertikalna struktura atmosfere iznad 
Beograda u danima sa kosavom 

Vertikalna struktura kosave je posebno proucavana tokom specijalnog 
osmatrackog perioda nazvanog ALPEX (Alpski eksperiment). Analiza 
prikupljenih podataka je pokazala da se maksimalni udari kosave javljaju skoro 
simultano sa formiranjem jezgra niske mlazne struje i da je pojava kosave 
pracena povecanjem stabilnosti atmosfere (Vukmirovic, 1985a). 

U ovom odeljku analizirana je vertikalna struktura troposfere u toku 
tridesetjednog dana neprekidnog duvanja vetra kosave (od 14. januara do 13. 
februara 1972. godine) izmerenog na stanici Beograd-Zeleno Brdo, zbog istih 
karakteristika udarnosti vetra kao u periodu od 1971. do 1975. godine. Ovo je 
bio rekordan broj dana sa kosavom u januaru i februaru od kada postoje 
meteoroloska osmatranja u Srbiji (od 1887. godine). Inace, srednji broj dana sa 
kosavom u januaru je 6,5, a u februaru 7,3 (Milosavljevic, 1950). 


5.1 Vertikalni profili vetra, temperature i 
relativne vlaznosti u danima sa kosavom 

Vertikalna struktura troposfere je predstavljena profilima srednjih vrednosti 
vetra, temperature i relativne vlaznosti do visine od 3 km u danima sa kosavom 
u januaru i februaru 1972. godine izmerenih na stanici Beograd-Zeleno Brdo, 
zajedno sa profilima srednjih desetogodisnjih vrednosti navedenih velicina za 
januar i februar racunatih u periodu od 1961. do 1970. godine (Unkasevic i 
saradnici, 1999). 
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Profili brzine vetra tokom januara i februara 1972. godine ukazuju na rast 
brzine sa maksimumom u sloju atmosfere izmedu 400 i 600 m iznad tla (slike 
25a i 26a). Iznad nivoa maksimuma vetra, centra niske mlazne struje, brzina 
vetra je brze opadala nego sto je prethodno rasla sve do visine od 1750 m iznad 
tla (850 mb), a zatim je ponovo bila u porastu. 

Ovakav profil brzine vetra pokazuje da se vetar kosava moze posmatrati 
kao specificna mezostruja, odnosno kao mlazna struja granicnog sloja atmosfere 
(Vukmirovic i Merkle, 1991). 

U daljoj analizi nije uocena bitna razlika izmedu profila brzine vetra kosave 
izmerene tokom januara i februara, narocito do visine od 750 m iznad povrsine 
tla. 

Vetar kosava znatno remeti srednje strujanje iznad Beograda. Naime, 
srednje brzine vetra u razmatranim mesecima 1972. godine u sloju do 1000 m 
iznad tla bile su vece od srednjih desetogodisnjih vrednosti vetra, a od oko 1250 
m bile su manje (slike 25a i 26a). 

Profil srednje temperature u 01 cas (po lokalnom vremenu) tokom januara 
ukazuje na postojanje izotermije (sloj vazduha u kome se temperatura ne menja 
sa visinom) u sloju do 400 m visine i inverzije (sloj vazduha u kome temperatura 
raste sa visinom) do visine od oko 750 m iznad tla (slika 25b). Danju, (u 13 
casova po lokalnom vremenu) zbog dnevnog zagrevanja i jake turbulencije, 
temperatura opada sa visinom nadadijabatski (1,2°C/100 m) do visine od oko 
550 m iznad tla. U sloju visine izmedu 550 i 1000 m iznad tla postojala je slaba 
inverzija. Iznad ovog sloja inverzije, vertikalni gradijent temperature bio je 
zanemarljiv (0,27 °C/100 m) do visine od oko 1250 m iznad tla. U sloju visine 
izmedu 1250 i 2750 m iznad tla vertikalni gradijent temperature iznosio je 
0,52°C/100m. 

Tokom meseca februara u 01 cas (po lokalnom vremenu) postojala je 
inverzija temperature do visine od oko 750 m iznad tla (slika 25b). Iznad ove 
visine, temperatura je opadala ravnomerno sa visinom sa gradijentom od oko 
0,51°C/100 m. 

Slicno kao i u mesecu januaru, u februaru je temperatura u 13 casova (po 
lokalnom vremenu) opadala sa visinom do visine oko 450 m iznad tla, nije se 
menjala do visine od oko 550 m iznad tla, da bi potom rasla sa visinom do oko 
750 m sa temperaturnim gradijentom u iznosu od 0,9°C/100 m. Vertikalni 
gradijent temperature iznad visine od 750 m iznosio je 0,51°C/100 m, tj ., kao i 
nocu. Poredenje ovih profila za januar i februar sa vertikalnim profilima 
temperature dobijenim za srednje desetogodisnje vrednosti temperature pokazuje 


da je u toku dana sa kosavom donja troposfera znatno stabilnija i toplija. 
Izuzetak je samo sloj visine do 400 m iznad tla u kome je, oko podne, postojao 
nadadijabatski gradijent temperature. 

Takode je uoceno da se prizemne inverzije mogu javiti i pri brzini vetra od 
16 m/s (2. i 3. februara 1972. godine u 01 cas), dok pridignute inverzije redovno 
prate kosavu danju (u 13 casova). Prizemne inverzije se pri kosavi uglavnom 
javljaju nocu i to u zimskom periodu. One se razvijaju ne samo pod uticajem 
procesa izracivanja i slabljenja turbulentne razmene, vec i drugih faktora kao sto 
je advekcija toplijeg vazduha u jugoistocnom strujanju, kontaktno hladenje 
vazduha od hladnije zemljine povrsine i drugi procesi (Vukmirovic i Merkle, 
1991). Debljina pridignutih inverzija u oba termina osmatranja tokom januara 
iznosila je izmedu 100 i 200 m. Tokom februara, debljina ovih inverzija iznosila 
je izmedu 300 i 400 m u 01 cas i izmedu 100 i 200 m u 13 casova, kao i u 
januaru. Ove pridignute inverzije kod nefrontalnih kosava su posledica 
vertikalnog spustanja vazduha u anticiklonu ili poticu od frontalne povrsine kod 
frontalnih kosava. 

Intenzitet prizemnih inverzija u januaru moze biti izmedu 4,5 i 5,0°C (u 01 
cas), dok su u proseku slabe i iznose izmedu 0,0 i 1,0°C. Intenzitet pridignutih 
inverzija je najcesce izmedu 1,1 i 3,0°C. 

Analiza vertikalnih profila vetra i temperature pokazuje da se nivo 
maksimuma brzine vetra nalazi ili na vrhu prizemne inverzije ili na bazi 
pridignute inverzije, odnosno na visini izmedu 400 i 600 m iznad tla. Ova 
cinjenica se mozda moze objasniti nemogucnoscu trenutnog uskladivanja polja 
vetra i temperature. 

Usled jakog vetra i intenzivnog mesanja vazduha u donjoj troposferi, 
relativna vlaznost vazduha se neznatno menja sa visinom. Naime, u mesecu 
januaru relativna vlaznost vazduha polako opada sa visinom do 750 m, a iznad 
ove visine u sloju do 3000 m iznad tla sporo se menja sa visinom u oba termina 
osmatranja (slika 25c). 

Tokom februara u nocnom terminu osmatranja (slika 26c) relativna vlaznost 
brzo opada do visine od 500 m, dok tokom dana se ne menja sa visinom, a zatim 
se do visine od 3000 m iznad tla skoro ne menja. Relativna vlaznost vazduha je 
nocu veca za oko 10% nego danju u oba posmatrana meseca. 

U poredenju sa vertikalnim profilom relativne vlaznosti vazduha dobijenim 
na osnovu srednjih desetogodisnjih vrednosti relativne vlaznosti (1961-1970), 
relativna vlaznost vazduha u danima sa kosavom u januaru bila je veca u oba 
termina osmatranja (slika 25c), dok je u februaru relativna vlaznost bila manja. 
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Jedini izuzetak je relativna vlaznost nocu u danima sa kosavom iznad visine od 
2000 m (slika 26c). 



Slika 25. Vertikalni profili: a) brzine vetra (m/s), b) temperature (°C) i c) relativne vlaznosti 
vazduha (%) racunati kao desetogodisnje srednje vrednosti u periodu od 1961. do 1970. godine 
(oznaceni sa 1) i u toku dana sa kosavom u januaru 1972. godine (oznaceni sa 2) izmereni na 

stanici Beograd -Zeleno Brdo. 
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5.2 Faktor ventilacije kod vetra kosave 

Pored udara kosave, koji nanose velike stete privredi, ovaj vetar dovodi do 
dobre provetrenosti i smanjenja zagadenosti prizemnog sloja atmosfere, narocito 
u gradskim i industrijskim oblastima. 

Za pozitivan ekoloski uticaj kosave bitna je njena moc ventilacije koja je 
data ventilacionim faktorom. Prema Summersovom modelu prostiranja zagade- 
nosti u atmosferi (Oke, 1978), koncentracija stetnih materija u vazduhu je 
obrnuto proporcionalna ventilacionom faktoru V e (m 2 s _I ) koji je definisan kao 
proizvod srednje brzine vetra v (m/s) u sloju ventilacije visine h (m), tj., kao: 

V e =vh. 

Pod slojem ventilacije ovde se podrazumeva sloj mesanja ili granicni sloj 
atmosfere. Sloj ventilacije kod kosave odreden je njenom vertikalnom 
strukturom, tj., vertikalnom raspodelom (profilom) brzine vetra i temperature 
(slike 25 i 26). U slucaju jakih vetrova, kao sto je kosava, visina tog sloja je 
jednaka nivou maksimalnog vetra u njegovom vertikalnom profilu, sto je 
najcesce na visini izmedu 300 i 400 m iznad tla. Medutim, visina sloja mesanja 
se neprekidno menja i nezavisno od doba dana ili noci krece se od 300 do 600 pa 
i 1000 m iznad tla doprinoseci tako velikoj vrednosti faktora ventilacije. 

Kao sto je istaknuto na pocetku ovog poglavlja, za kosavu je tipicno da 
najcesce duva u sloju atmosfere ogranicenom pridignutom inverzijom 
temperature, cija se baza najcesce podudara sa maksimalnom brzinom u 
vertikalnom profilu vetra. To ja za kosavu prirodan i dobro definisan 
ventilacioni sloj, sto ima za posledicu veoma delotvoran faktor ventilacije. 

Treba naglasiti da pojava prizemnih inverzija temperature ne smanjuje moc 
ventilacije kosave, jer ne utice bitno na visinu sloja mesanja, a kao sto je vec 
istaknuto, one mogu da postoje i pri brzinama vetra od 16 m/s na visini 
anemografa. 

U tabeli 14 date su vrednosti ventilacionog faktora za slucaj kosave koja se 
javila u periodu od 18. do 24. marta 1982. godine u Beogradu-Opservatorija 
(Vukmirovic, 1997). Iz tabele 14 se vidi da faktor ventilacije, V e9 ima najvece 
vrednosti pri velikim srednjim brzinama kosave ili usled porasta visine sloja 
mesanja preko 400 m, koji je obicno kratkotrajan. Najvece brzine i udare pri 
zemlji kosava ima kada se na nivou maksimalne brzine javi tzv. jezgro mlaza. 
Pomenuta tabela takode pokazuje da ventilacioni faktor ima visoke vrednosti 
kako danju tako i nocu, bez obzira na brzinu kosave pri zemlji v (m/s). 


Tabela 14. Vrednosti faktora ventilacije za kosavu izmerene na stanici Beograd-Opservatorija 
(18.-24. mart 1982. godine, Vukmirovic 1997). 


Dan 

Cas 

V 

K 

Dan 

Cas 

V 

K 

Dan 

Cas 

V 

v. 

18. 

12:15 

2 

4020 

21. 

06:15 

10 

5850 

22. 

21:40 

7 

9900 


18:15 

9 

8400 


09:35 

12 

3500 

23. 

00:50 

7 

12180 

19. 

00:15 

14 

7650 


12:15 

11 

3300 


03:15 

11 

12470 


10:17 

12 

5600 


15:15 

5 

3525 


07:25 

9 

20020 


18:19 

6 

5400 


18:15 

7 

4600 


09:15 

7 

20375 

20. 

00:15 

9 

12400 

22. 

00:25 

7 

6625 


12:15 

10 

21245 


06:30 

4 

14500 


03:15 

6 

4500 


15:35 

5 

13920 


12:30 

J 

9450 


06:15 

7 

6600 


18:19 

6 

16500 


18:15 

8 

11700 


09:15 

8 

4500 

24. 

00:15 

8 

15000 


22:00 

5 

9360 


12:30 

8 

4950 


06:37 

4 

5060 

21. 

00:15 

5 

9100 


15:15 

8 

4500 


12:34 

3 

2700 


03:15 

5 

5500 


18:15 

7 

9375 


19:00 

4 

950 


Ako se ove vrednosti faktora ventilacije uporede sa faktorom ventilacije 
racunatim za pretpostavljene srednje vrednosti visine mesanja i srednje brzine 

vetra u njemu (oko 200 m 2 s _l za noc i oko 1500 m 2 s _l za dan) vidi se koliko je 
znatna moc provetravanja pri pojavi kosave. 

O prostornoj raspodeli faktora ventilacije moze se suditi samo 
kvantitativno, na primer, iz anemografskih zapisa udarnosti u danima sa 
kosavom. Jedna takva analiza epizode jake kosave Popovica (1982) pokazuje da 
je najveci ventilacioni efekat kosave u istocnom delu kosavske oblasti tj., u 
zavetrini Karpatskog planinskog sistema (slika 27). 
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Slika 27. Izotahe (izolinije maksimalnih udara vetra) kosave (km/h) u periodu od 16. do 18. 
februara 1979. godine (Popovic, 1982). 


6 Spektralna analiza vetra kosave 

Spektar vremenske serije (VS) predstavlja raspodelu varijanse serije kao 
funkciju ucestanosti. Ocena i proucavanje spektra je cilj spektralne analize. 
Postavlja se pitanje zasto analizirati spektar? Spektar je od interesa jer mnoge 
prirodne pojave zavise od ucestanosti, a razumevanje zavisnosti od ucestanosti 
moze da nas informise o fizickim mehanizmima, koji bitno odreduju posmatrani 
meteoroloski niz. 

Do znacajnog saznanja o prirodi vremenskih serija mozemo doci ako 
izvrsimo transform aciju autokorelacione funkcije (ACF). Ova operacija 
transformacije autokorelacione funkcije iz domena vremena u domen ucestanosti 
daje novu funkciju koju zovemo spektar varijacija. Spektar i ACF su dva 
komplementarna nacina prikazivanja varijanse VS, jer pruzaju istu informaciju 
ali izrazenu na razlicite nacine. Autokorelaciona funkcija nam daje informacije u 
wemenskom domenu , a spektar u domenu ucestanosti ili frekvencij shorn domenu. 
U nekim situacijama (istrazivanje vulkanskih erupcija) je korisniji pristup u 
vremenskom domenu zasnovan na ACF, dok je u drugim situacijama (na primer, 
sezonske promene povrsinske temperature) bolji pristup u domenu ucestanosti 
zasnovan na spektru. 

Neparametarske ili klasicne metode ocene spektra razvijene su u cilju 
poboljsanja ocene spektra pomocu periodograma. Poboljsanja se ogledaju u 
tome sto se na periodogram ili na autokorelacionu funkciju primenjuju operacije 
osrednjavanja, tako da se dobije konzistentna ocena spektralne gustine. Pri tome 
se ne cine nikakve pretpostavke u pogledu karakteristika posmatranog niza. 

U ovoj Glavi prikazace se Blekman-Takijev metod, metod sa vise prozora i 
vejvlet transformacija (Tosic, 2005). 
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6.1 Matematicki opis spektra primenom 
Blekman-Takijevog metoda 


Blekman-Takijev (BT) metod (Blackman i Tukey, 1958; Tosic, 2004) je 
jedan od neparametarskih metoda za ocenu spektra. 

Ako imamo N clanova serije , izracunamo serijske kovarijanse C T za 
sve pomake odx = 0 do x - m (m < N) prema: 

| /V-T 

C T = — — - Yj ( x i - X )( X i + r - x ) > 

N " 1 /=1 


gde je x srednja vrednost od svih x t . Spektralne ocene S k , k = 
se zatim dobijaju pomocu jednacina: 

0,1,2,..., m-\,m 

1 1 W_1 

*^o ~ 7 (Q + 9 

2m m T=1 

(1) 

- C 0 2 ^ ( nki\ 1 r , n k 

S k - — + — >C x cos + C m (-1) , 

m m \ m ) m 

(2) 

i r 11 W_1 

(3) 


Jednacina (1) se koristi za izracunavanje spektralne ocene S 0 , koja 
odgovara beskonacnoj talasnoj duzini (trendu), a jednacina (3) se koristi za 
izracunavanje spektralne ocene S m , koja odgovara najkracoj talasnoj duzini 
resivoj u spektru (talasnoj duzini jednakoj dvostrukom intervalu izmedu 
uzastopnih osmatranja). Preostalih m - 1 vrednosti spektralnih ocena se 
izracunavaju iz jednacine (2) stavljajuci da k u kosinusnom argumentu uzima 
uzastopne vrednosti k = 1, 2,..., m - 1 . 

Konacne spektralne ocene, S k9 se zatim izracunavaju uravnavanjem 
spektralnih ocena pomocu formula: 

S 0 =i(s 0 +<$,), (4) 

S/c =“ (^*-1 + 2 k< + Sk + 1 )> 


(5) 
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S m ~ ~ (^m-1 +£/;;)> ( 6 ) 

gde k u jednacini (5) uzima vrednosti 1, 2, 1. To je tzv. Hamingov 

metod. 

6.1.1 Test znacajnosti koji se primenjuje na spektre 

Znacajnost pikova spektra se ispituje u odnosu na tzv. nulti spektar, koji 
moze da bude “beli sum” ili “crveni sum”. Ako se koeficijent serijske korelacije 
nakon prvog koraka r, ne razlikuje od nule za statisticki znacajan iznos, nulti 
spektar je “beli sum”, ili drugim recima horizontalna prava linija cija je vrednost 
svuda jednaka srednjoj vrednosti od svih m + 1 neuravnatih spektralnih ocena u 
izracunatom spektru. 

Medutim, ako se koeficijent serijske korelacije razlikuje od nule za 
statisticki znacajan iznos, treba proveriti da li koeficijenti korelacije priblizno 
prate relaciju 

r 2 « r, 2 ; r 3 » r , 3 ; itd. 

U tom slucaju, odgovarajuca nulta hipoteza treba da bude “Markovljev 
crveni sum”. 

Prema Gilmanu i sar. (1963), spektar “crvenog suma” je: 


P k =s 


+ r? - 2 r { cos(7r&/m) 


gde je s srednja vrednost od svih m + 1 spektralnih ocena S k u izracunatom 
spektru. Kada je r, « 0 u (7), dobija se spektar “belog suma”. 

Statistika koja je povezana sa svakom spektralnom ocenom, na osnovu koje 
se odreduje znacajnost odstupanja te ocene od nultog spektra, je odnos uzoracke 
spektralne ocene ( S k ) i nultog spektra. Tukey je (1950) otkrio da je ova velicina 
raspodeljena kao hi-kvadrat podeljena sa brojem stepeni slobode (v), koji su 
bazirani na duzini serije od N vrednosti i na maksimalnom pomaku od m 
vremenskih jedinica: 


Kriticni nivoi % 2 / v raspodele za odredeno v se dobijaju iz tabela za % 2 . 
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6.2 Metod sa vise prozora 

Thomson je 1982. godine uveo metod sa vise prozora (MTM, engl. 
Multitaper method). Kod ovog metoda su podaci pomnozeni ne sa jednim, vec 
sa nekoliko prozora koji su otporni na spektralno curenje. Ovo rezultuje u 
nekoliko ( S ) taperovanih vremenskih serija, a s (t)x(t), od jedne pocetne 
vremenske serije x (t) , gde s = 1 Prozori su ortogonalni, tako da svaki od 
njih uzorkuje vremensku seriju na drugaciji nacin. Statisticka informacija 
odbacena prvim prozorom je delimicno pokrivena drugim prozorom, informacija 
odbacena sa prva dva prozora je delimicno pokrivena trecim prozorom, itd. 
Primenjuje se samo nekoliko prozora nizeg reda jer spektralno curenje raste sa 
porastom reda. Primenjujuci Furijeovu transform aciju na svaku od S 
taperovanih vremenskih serija, proizvodi se S ocena spektra, Y s (f) . Prozori 
a s se nazivaju “sopstveni prozori” jer nastaju razlaganjem na sopstvene 
vrednosti. Slicno se i spektri Y s (/) zovu sopstveni spektri. Spektar dobijen 
koriscenjem vise prozora ima manju varijansu od jednom taperovanog spektra i 
u isto vreme je otporan na spektralno curenje. 

Za datu vremensku seriju x(t ), odredujemo prvo set od S ortogonalnih 
Slepian prozora a s (t) i njihovih S taperovanih Furijeovih transformacija ili 
sopstvenih spektara Y s (/) , s = 1,..., S , kao: 

n(/) = Z^(040e ,2,t/ ' A ', (9) 

t = \ 

gde je a s (t) 5-ti clan u familiji ortogonalnih Slepian prozora, sa 5 = 1,..., S, a 
At je interval uzorkovanja. Mogu se resiti samo spektralne fluktuacije na 
ucestanostima vecim od Rejlijeve frekvencije. Frekvencijska rezolucija spektra 
je data sa 2 pf R , gde je p (ceo broj) parametar vremensko-frekvencijskog 
propusnog opsega (engl. bandwidth). Prema Park i sar. (1987), samo prvih 
S = 2p-\ prozora je otporno na spektralno curenje. U klimatskim 
proucavanjima, izbor p = 2 i stoga S = 3 obezbeduje dobar kompromis izmedu 
rezolucije koja odgovara resavanju klimatskih signala i varijanse spektralne 
ocene (Mann i Park, 1994, 1996; Tosic i Unkasevic, 2004; Tosic i Unkasevic, 
2005). 


Linearna kombinacija sopstvenih spektara data sa Y (/) obezbeduje spektar 
sa optimalnim osobinama kao kompromisno resenje izmedu spektralne 
rezolucije i varijanse: 

ixiw)i 2 

^(/) =i =S — • (>°) 

5=1 

Znacajnost periodicnih ili kvazi-periodicnih signala se meri u odnosu na 
pozadinski sum. Nivoi znacajnosti ocene suma se izracunavaju pretpostavljajuci 
da spektar ima xl raspodelu sa v«2 S stepeni slobode (Mann i Lees, 1996). 
Pretpostavlja se da je odnos snage pika i lokalnog nivoa suma u spektru 
raspodeljen kao xV v 1 moze se porediti sa tablicnim % 2 vrednostima 
verovatnoce raspodele da bi se odredila znacajnost pika. 


6.3 Vejvlet transformacija (WT) 


Vejvlet analiza (WA) je metod za analizu promena snage VS sa vremenom. 
Vejlet transformacijom se vremenska serija transformise iz vremenskog u 
vremensko-frekvencijski prostor. Nasuprot Furijeovoj transformaciji, koja 
omogucava prikaz globalnog frekvencijskog sadrzaja, WT obezbeduje 
lokalizovanu, trenutnu ocenu amplitude i faze za svaku spektralnu komponentu 
serije. Ovo daje prednost WA u analizi nestacionarnih podataka kod kojih 
amplitude i faze harmonijskih komponenti mogu da se brzo menjaju u vremenu 
ili prostoru (Torrence i Compo, 1998). 

Vejvleti su familije funkcija koje se generisu skaliranjem i translacijom 
jedne prototipske funkcije (engl. mother wavelet) G(t) : 


G a +(t) = -=G 
da 


t-b 
a ) 


( 11 ) 


Prototipska vejvlet funkcija mora da zadovolji sledeci uslov 

\G(t)dt = 0 , 


( 12 ) 


sto ukazuje na njen oscilatorni karakter, zbog cega je ova funkcija i dobila ime 
talasic (Mojsilovic, 1997). Realan skalar b se naziva parametar translacije i 
odgovara centralnoj tacki vejvleta u VS. Realan i pozitivan skalar a je 
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parametar skaliranja i odreduje sirinu vejvleta. Faktor 1 /yfa normalizuje 
vejvlete tako da imaju jedinicnu energiju i otuda su uporedivi za sve razmere a . 

Osnovu WT cini ideja da se proizvoljna funkcija x(t) predstavi kao 
superpozicija vejvleta. Drugim recima, WT predstavlja dekompoziciju funkcije 
x na sledeci nacin: 

1 

Ukoliko se vrsi dekompozicija realnih signala, znak za konjugovano 
kompleksnu vrednost * se moze izostaviti. Prethodna jednacina predstavlja 
definicioni izraz za kontinualnu vejvlet transformaciju. Jednacina (13) opisuje 
postupak dekompozicije prilikom izracunavanja WT (formula za analizu). 
Sinteza funkcije x(t) , na osnovu transformacionih (vejvlet) koeficijenata 
T(a,b) dataje sledecim izrazom: 


T(a,b) = (x(t ), G a b (/)} 


d/,aeR + ,6eR. (13) 


x(t) = c j \ T(a,b)G a<b (t)^- 


gde se pretpostavlja da su i x(t) i G ab (t) realne funkcije. Da bi vazila 
rekonstrukciona formula (14), vejvlet G ab (t) mora da zadovolji uslov (12). 


Jedan od najcesce koriscenih vejvleta u atmosferskim i geofizickim 
naukama je Morletov vejvlet: 

/ 2 

G(t) = e ^e lc ' , (15) 


koji se sastoji od ravnog talasa e' c/ frekvencije c-f (ili talasnog broja c = k u 
prostornom domenu), modulisanom u vremenu Gausovom envelopom jedinicne 

sirine e ' . 

Praktican postupak za izvrsenje vejvlet analize se moze sumirati na sledeci 
nacin. Kako radimo sa vremenskim serijama konacne duzine, dogadaju se greske 
na pocetku i na kraju vejvlet spektra, jer Furijeova transformacija (FT) 
pretpostavlja da su podaci ciklicni. Zbog toga treba taperovati vremensku seriju 
x(t) koristeci prozor oblika 1-sincp, koji tezi nuli ka krajnjim tackama serije. 


Konacan broj tacaka u seriji treba da bude stepen od 2, tj. N = 2 m , gde je m ceo 
broj da bi lakse izracunali brzu Furijeovu transformaciju (FFT). Zatim se 
oduzme srednja vrednost taperovane vremenske serije i izracuna se njena 
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Furijeova transformacija x(f) . Takode, treba izracunati i FT vejvleta G(t) za 

specifican razmer a da bi dobili G(af). Moze se koristiti proizvoljan broj 
razmera u obliku: 

a t = a 0 2 J&J , j = 

gdeje 


7 = (8,riog/^\ 

V a 0 ) 

Najmanji razmer a 0 je odreden sa a 0 = 2 St , a najveci razmer je odreden sa 
J (Torrence i Compo, 1998). Vremenski interval je odreden sa 5/ , a manje 
vrednosti 5 j daju bolju rezoluciju. U nasem slucaju izabrano je 5 j = 0,125 . 

Zatim se izracunava WT kao u jednacini (13) konvolucijom proizvoda 
G (af)x(f) u Furijeovom prostoru. Konacno, uradi se inverzna FT da bi dobili 
T(b,a) kao funkciju vremena b i razmera a. Posto je vejvlet funkcija G(t) 
kompleksna, rezultujuca vejvlet transformacija je takode kompleksna. Stoga se 
moze napisati kao suma realnog Re{r(Z?,^)} i imaginarnog Im{r(Z?,^)} dela, ili 
kao kombinacija amplitude i faze. Dodatno, mozemo definisati “vejvlet spektar 
snage” kao |r(Z?,a)| 2 . WT mozemo predstaviti crtanjem vejvlet spektra kao 
funkcije b na linearnoj x-osi i razmera a na logaritamskoj j;-osi ili u funkciji 
Furijeove periode i vremena (tj. parametra translacije). Medutim, ponekad se 
moze naci veza izmedu talasne duzine i razmera a. Furijeova perioda P je 
povezana sa parametrom skaliranja a i parametrom Morletove vejvlet funkcije 
co 0 pomocu sledece relacije (Torrence i Compo, 1998): 



Dopunjavanje nulama uvodi diskontinuitete u krajnjim tackama, pa kada se 
ide ka vecim razmerama, opada amplituda blizu ivica. Zbog toga se uvodi oblast 
uticaja konusa (COI), a bira se tako da vejvlet spektar opada za faktor e~ 2 , sto 


obezbeduje da su uticaji ivica zanemarljivi izvan ovih tacaka. 
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6.4 Analiza rezultata 

U ovom odeljku su prikazani rezultati dobijeni koriscenjem Blekman- 
Takijevog (BT) metoda, metoda sa vise prozora (MTM) i vejvlet transformacije 
za kosavu sa najduzim trajanjem (14. januar-13. februar 1972. godine). Za prva 
dva metoda na x-osi je prikazan period u casovima, a na_y-osi spektralna gustina. 
Za treci metod je na jc-osi prikazano vreme u casovima, na ^-osi period u 
casovima, a treca promenljiva je vejvlet spektar. 

Ono sto je odmah uocljivo je maksimum sa periodom od oko 120 h 
znacajan na 95% nivou, dobijen pomocu sva tri metoda. Sa slika 28 i 29 se vidi 
da su spektri dobijeni BT i MT metodom veoma slicni. 

Koristeci Blekman-Takijev metod (slika 28a), dobijene su oscilacije od 
121,3 h, ali i krace od 4,0; 3,6; 3, 3-3, 2; 3, 0-2, 9; 2, 5-2, 3 h za srednje casovne 
vrednosti, dok su za maksimalne casovne vrednosti kosave dobijene oscilacije 
od 121,3; 7, 6-7, 5; 5,7; 5,4; 4,8; 4, 4-4, 3; 3,0; 2,7; 2,5 i 2,3 h (slika 28b). 

Koristeci MTM (slika 29a), dobijene su oscilacije od 102,4-128; 4,4; 3,0; 
2,9; 2,3 za casovne, a za maksimalne casovne vrednosti kosave su dobijene 
oscilacije od: 102,4; 4,0; 3,6; 2,9, 2, 3-2, 5 h (slika 29b). 

Koristeci vejvlet transformaciju, dobija se organizovana struktura od oko 
120 h i za casovne i za maksimalne satne vrednosti vetra, ali je za maksimalne 
casovne vrednosti 95% kriva uza (slika 30). Ta perioda od oko 120 h se javlja u 
prvih 400 casova. Takode se uocavaju i manje oscilacije od 2 do 7 h i to u toku 
prvih 200 casova i od 400-tog do 600-tog casa. 

Iz prethodne analize se moze zakljuciti da je dominatan period od oko 120 
h, tj., oko 5 dana, sto nam ukazuje da je kosava pojava sinoptickih razmera, sto 
je vec uoceno na osnovu radiosondaznih podataka (ALPEX). 
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k) Period (sati) 

Slika 28. Blekman-Takijev metod za: a) srednje casovne i b) maksimalne casovne vrednosti 
kosave (14. januar-13. februar 1972. godine) registrovane na Opservatoriji u Beogradu. 
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7 Modeliranje vetra kosave 

U ovoj Glavi monografije prikazace se modeliranje dve sinopticke situacije 
u kojima je duvala kosava: od 05. do 06. decembra 1995. godine i od 14. do 17. 
aprila 1994. godine. Prvu sinopticku situaciju, razmatranu u Glavi 3, je 
karakterisala jaka kosava. Druga sinopticka situacija je pracena slabom kosavom 
i transportom saharske prasine. 

7.1 Modeliranje sinopticke situacije iznad 

Evrope u periodu od 5. do 6. decembra 1995. 
godine 

Za simulaciju kosave Lazic i Tosic (2000) su koristili model za ogranicenu 
oblast sa stepenastim planinama (Mesinger, 1984), tzv. Eta model. Model je 
definisan na polurazmaknutoj Arakavinoj E mrezi. Ima ugradenu specijalnu 
tehniku za sprecavanje razdvajanja resenja (Mesinger, 1973; Janjic, 1979). 
Vertikalna kooordinata je eta koordinata, kao general izacij a sigma koordinate 
(Philips, 1957), sa stepenastom prezentacijom planina. U horizontalnoj advekciji 
model ima ugradenu kontrolu nelinearne kaskade energije (Janjic, 1984). 
Vremensko diferenciranje je eksplicitno sa rascepljivanjem (Mesinger, 1977; 
Janjic, 1979). 

Fizicki paket modela ukljucuje Melor-Jamada parametrizacionu semu nivoa 
2,5 za planetarni granicni sloj i Melor-Jamada parametrizacionu semu nivoa 2 za 
prizemni sloj. Parametrizovani su: povrsinski procesi, turbulencija, padavine 
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velikih razmera, konvekcija, horizontalna difuzija i zracenje (Janjic, 1990; 
1994). 

Za eksperimente je koriscen model sa horizontalnom rezolucijom 
0,25° x 0,25° (28 km x 28 km) i 16 slojeva po vertikali. Elementarni vremenski 
korak je bio 90 s. Horizontalna oblast modela je bila definasana izmedu 0° i 30° 
E geografske duzine i od 40° do 50° N geografske sirine (slika 31). 

Posto je Eta model, model za ogranicenu oblast, potrebno je da se u prvom 
spoljnom redu tacaka koriste vrednosti dobijene integracijom nekog modela cija 
je integraciona oblast veca od oblasti modela koji se posmatra. Stoga su 
vremenski zavisne granicne vrednosti uzete iz analiza Evropskog centra za 
srednjorocnu prognozu vremena (ECMWF), koje su linearno interpolisane iz 
analiziranih polja raspolozivih na 6 sati. Kako integraciona oblast modela 
pocinje od treceg reda tacaka, vrednosti u drugom redu se dobijaju 
kombinacijom vrednosti iz prvog i treceg reda. 

Tipicna kosavska situacija se dogodila od 5-7. decembra 1995. godine (slika 
31) i detaljno je opisana u Glavi 3. Koristeci simuliran vetar, izracunate su 
trajektorije unapred polazeci od osam najnizih nivoa modela (LT=1 sa 
pribliznom visinom od z=3820 m; LT=2, z= 3057 m; LT=3, z=2380 m; LT=4, 
z=1784 m; LT=5, z=1264 m; LT=6, z=818 m; LT=7, z=442 m i LT=8, z=136 m) 
za izabranu pocetnu poziciju A,=24°, cp=44° i pokazane su na slici 32. 

Trajektorija koja je krenula sa najnizeg nivoa modela (LT=8 na z=136 m ili 
ri=0,98) je bila blokirana planinanama (stepenasta eta orografija od 272 m) u 
pocetnom polozaju i nije prikazana na crtezu. Sve trajektorije povecavaju brzinu 
u pocetku integracije. Maksimalna brzinaje bila oko 3 h od pocetka integracije, 
u orografiji tipa kanala Iron Gate kroz Karpate. Sve trajektorije su se penjale u 
pocetku iznad orografskih stepenika od 613 m, 1023 m, a neke od njih i iznad 
1506 m, a zatim su se spustale u zavetrenoj strani Karpata. Mlazna struja je 
simulirana na oko 1200 m sa brzinom vetra od 25 m/s. U isto vreme 
maksimalna brzinaje bila 21 m/s na nivou od 800 m, a maksimalna brzinaje 
bila 22 m/s na nivou od 1 800 m. 

Trajektorije koje su krenule sa nivoa 1 i 2 (LT=1, z=3820 m i LT=2, z=3057 
m) su imale skoro jugoistocni pravac, povecavajuci visinu iznad stepenica od 
1506 m preko Karpata (slika 32 a, b). Trajektorija koja je krenula sa nivoa 3 
(LT=3, z=2380 m) je dostigla maksimalnu brzinu od 18 m/s nakon 3 h od 
integracije, iznad orografske stepenice od 1506 m (slika 32 a, c). Trajektorija 
koja je krenula sa nivoa 4 (LT=4, z=1784 m) je imala maksimalnu brzinu od 22 
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m/s nakon 3 h od pocetka integracije, upravo iznad orografskog stepenika od 
1506 m. 



O -WOE i-00 0.1M7E+02 

Minimum Vector Mnximum Vector 



0 2R-JE+00 


0.23GE+02 
Maximum Vector 

Slika 31. Simulirana vremenska situacija za 05. decembar 1995. godine na 850 hPa: a) 00 h. 


b) +24 h, c) +48 h. 
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Trajektorija koje je krenula sa nivoa 5 (LT=5, z=1264 m) je imala 
maksimalnu brzinu od 25 m/s u trecem satu integracije. To je niska mlazna 
struja vetra kosave, locirana oko 1200 m u vertikalnom pravcu, sa jugoistocnim 
pravcem, u kanalu kroz Karpate. Trajektorija koja je krenula sa nivoa 6 (LT=6, 
z=818 m) dostigla je maksimalnu brzinu od 21 m/s nakon 3 h. Trajektorija koja 
je krenula sa nivoa 7 (LT=7, z=442 m) je imala maksimalnu brzinu od 12 m/s 
nakon 4 h. Tipicna raspodela brzina vetra duz trajektorija se poklapa sa 
osmatranjima. Svim trajektorijama opada brzina vetra, nakon prolaska kosavske 
oblasti. 


7.2 Modeliranje sinopticke situacije iznad 
Evrope u periodu od 14. do 17. aprila 1994. 
godine 


0.3t 

0.4; 

0.5 

O.67 

I 

0.7- 

0 . 8 : 

0.9; 

1.0 


^ # o' 

c .o 


r ^ <> 

, c -' A * D ‘ 



4 o n 

a A O a 

a B “ 

C>o o ^ o xr 

. ° ° ° ■ * T 


° r — 


(b) 


0 3 6 9 12 15 18 21 24 27 30 33 36 

t(h) 



Slika 32. Trajektorije unapred za 05. decembar 1995. godine: a) horizontalni polozaj, 
b) vertikalan polozaj, c) brzina vetra duz trajektorija. 


Primer umerene kosave u periodu od 14-17. aprila 1994. godine je opisan u 
Vukmirovic i saradnici (2004). Balkansko poluostrvo je bilo pod uticajem 
transporta i depozicije saharske prasine. Ta prasina je bila bogata gvozdem koje 
daje crvenu ili zutu boju padavinama (Vukmirovic i saradnici, 1997). Razvoj 
depresije u severnoj Africi je poceo 14. aprila 1994. godine. Sinopticka situacija 
od 14. do 17. aprila je simulirana koristeci Eta model sa horizontalnim 
rezolucijama od 1 x 1° i 1 0 f x 10*. Slika 33a, b pokazuje 72 h prognozu vetra i 
geopotencijalne visine na 850 hPa koristeci manju i vecu rezoluciju, respektivno. 
Prasina podignuta u ciklonskoj cirkulaciji, nosena je ka jugoistocnoj Evropi i ka 
Beogradu. Padavine su se prvo javile u juznoj Srbiji tokom noci 15/16 aprila. 
Zuta kisa je zabelezena 16. i 17. aprila, kada su u Beogradu bile maksimalne 
padavine (Vukmirovic i saradnici, 1997). Kosava je bila umerene jacine tokom 
ovog slucaja saharske prasine. Maksimalna brzina vetra od 4,6 m/s je izmerena u 
24 h intervalu od 16. do 17. aprila. Nakon tog intervala vetar je bio promenljivog 
pravca sa brzinama ispod 1 m/s. 

Trajektorije unazad, izracunate koristeci prognozu iz Eta modela za 72 h sa 
veoma visokom rezolucijom, su prikazane na slici 34. Trajektorije na vecim 
visinama su brzo napustile oblast integracije, koja je bila mala u poredenju sa 
domenom integracije za rezoluciju od l°xl 0 . Slika 34b pokazuje vertikalnu 
poziciju trajektorija prema eta nivoima. Prilazeci Beogradu, trajektorije sa nizih 
nivoa imaju SE pravac (slika 34a) sa malim brzinama vetra (slika 34c). 
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Slika 33. 72h prognoza vetra i geopotencijalne visine na 850 hPa za 14 april 1994. godine 


koristeci rezolucije: a) l°x 1°, b) 10' x 10'. 
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0 LT=1 O LT=2 * LT=3 $ LT=4 □ LT=5 A LT=6 
(z=5 129 m) (z=2063 m) (z=1423 m) (z=707 m) (z=545 m) (z=347 m) 


eta 



t(h) 

Slika 34. Trajektorije unazad za 14. april 1994. godine: a) horizontalni polozaj (brojevi duz 
trajektorija oznacavaju vreme integracije), b) vertikalan polozaj, c) brzina vetra duz trajektorija. 
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Umerenom kosavom je transportovan vazduh kroz basen reke Morave 
izbegavajuci jako zagaden rudnik u Boru i basene visokih peci u istocnoj Srbiji. 
Dominantni metali iz Bora su Cu i Zn. Usled umerenih brzina vetra nije bilo 
uslova za resuspenziju delica iz dubokih dolina u atmosferu (Braaten et al, 
1990). Treba napomenuti da glavni izvori metala iz zelezara u centralnoj Srbiji 
nisu radili u to vreme. Prema analizi tri najnize trajektorije, transport metala iz 
Makedonije i juzne Srbije kosavom je bio dominantan u osmotrenoj epizodi. 
Pretpostavlja se da su frakcije Pb i Cd, nadene u uzorcima suve depozicije, 
transportovane iz severozapadne Makedonije, preciznije iz industry skog grada 
Tetova. Turbulentni tok je pojacao proces koagulacije u pocetku cistih delica 
prasine sa delicima koji sadrze Pb i Cd. Prilazeci Beogradu sve trajektorije se 
poklapaju i proces koagulacije je verovatan u poslednjih 6 h (slika 34). Izgleda 
da se ovaj proces nastavio i u severnoj Srbiji proizvodeci kontaminirano plodno 
zemljiste sa toksicnim metalima ili depozicijom ili resuspenzijom. 
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